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ABSTRACT 


This  study  is  concerned  with  the  geocentric  attitude  stabilization  of  an  earth 
satellite.  The  system  considered  in  this  report  utilizes  differential  gravity  to 
control  the  vehicled  orientation  and  gyroscopes  to  provide  vehicle  damping.  Several 
different  gyroscope  configurations  that  provide  indirect  gyroscope  damping  about 
all  three  of  the  vehicle  axes  are  presented.  For  each  of  these  configurations,  the 
equations  governing  the  dynamics  of  the  system  are  systematically  developed  in  terms 
of  the  attitude  deviation  angles  and  the  parameters  of  the  orbit.  The  expressions 
for  the  steady  state  deviation  angles  and  the  torque  which  must  be  supplied  by  gyro¬ 
scope  torque rs  are  also  obtained. 

The  analytical  procedure  utilized, considers  the  satellite  vehicle  to  be  composed 
of  a  rigid  body  to  which  is  affixed  internal  gyroscopes.  The  dynamics  of  an  orbiting 
rigid  body  are  first  obtained  with  the  inclusion  of  the  differential  gravity  torque. 

This  is  followed  by  a  development  of  the  equations  of  dynamics  of  a  gyroscope  con¬ 
strained  in  an  arbitrary  position  within  an  orbiting  body.  The  equations  obtained  for 
the  gyroscope  and  the  rigid  body  are  completely  general  in  that  they  are  valid  for 
an  elliptical  orbit.  The  resulting  equations  of  a  gyroscope  and  a  rigid  body  are 
combined  to  obtain  both  the  statics  and  dynamics  of  an  orbiting  vehicle  containing 
internally  mounted  gyroscopes.  This  analysis  is  restricted  to  a  circular  orbit.  To 
obtain  a  better  understanding  of  the  mechanism  of  gyroscope  damping,  an  analysis 
of  a  simple  configuration  employing  a  gyroscope  for  damping  purposes  is  also 
included. 
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angular  velocity  of  the  gyroscope  with  respect  to  inertial  space. 

set  of  Eulerian  angles  relating  the  reference  coordinate  system  to 
the  body  coordinate  system. 

set  of  Eulerian  angles  relating  the  body  coordinate  system  to  the 
gyroscope  coordinate  system. 

first  and  second  derivatives  with  respect  to  time. 


Page 

Appendix  A 

Derivation  of  Euler*  s  Equations 

106 

Appendix  B 

Derivation  of  Gyroscope  Equations 

110 

Appendix  C 

Derivation  of  Gravity  Gradient  Torque 

114 

Reference 

1 


Chapter  I  -  INTRODUCTION 

There  are  numerous  applications  which  require  a  particular  axis  of  a 
satellite  vehicle  to  be  oriented  in  some  specific  direction.  The  need  for  space 
vehicle  attitude  control  and  stabilization  can  only  be  appreciated  when  consideration 
is  given  to  the  mission  that  the  vehicle  is  to  perform.  For  example,  conversion  of 
solar  energy  to  electrical  energy  will  require  solar  energy  gathering  devices  which 
must,  by  necessity,  be  oriented  toward  the  sun.  Midcourse  maneuvers  requiring 
the  firing  of  rocket  engines  would  also  require  vehicle  attitude  control  to  give  the 
proper  direction  to  the  thrust.  As  the  satellite  moves  in  orbit  around  the  earth  it  might 
be  desired  to  obtain  weather  information  or  obtain  surveillance  information  over 
enemy  territory.  The  latter  applications  require  a  camera  to  be  so  oriented  that 
its  axis  is  aligned  with  the  geocentric  vertical  of  the  earth. 

The  various  sources  of  attitude  disturbance  torques  for  a  vehicle  orbiting  the 
earth  are  worthy  of  consideration.  A  brief  discussion  of  disturbance  torques  is  given 
below: 

a)  Aerodynamic  Torque 

At  an  altitude  of  about  500  miles  a  very  sparce  atmosphere  exists. 
Nevertheless,  the  torque  produced  by  aerodynamic  drag  can  be  significant.  This 
would  be  particularly  true  if  the  space  vehicle  was  to  employ  large  solar  energy 
collecting  panels  where  the  location  of  the  center  of  aerodynamic  pressure  might  not 
coincide  with  the  center  of  mass  of  the  vehicle-  The  torque  produced  is  a  function  of 
the  angle  of  attack  of  the  vehicle  and  for  a  constant  angle  of  attack  and  a  circular 
orbit,  the  torque  produced  will  be  constant. 

At  an  altitude  of  about  500  miles,  the  dynamic  pressure  is  approximately 
6xl0“5  dyne/ cm2.  For  a  space  vehicle  employing  flat  solar  energy  collecting  panels, 
100  feet  2  in  area,  with  a  2.  5  inch  offset  of  the  center  of  aerodynamic  pressure 
from  the  mass  center  of  the  vehicle,  a  disturbance  torque  of  approximately 
56  dyne-cm.  will  result  for  the  worst  angle  of  attack.* 

b)  Solar  Radiation  Pressure 

Solar  radiation  pressure  is  caused  by  the  impact  of  photons  with  the 
surface  of  the  vehicle.  As  in  the  case  of  aerodynamic  drag,  torques  produced  can  be 
significant  if  large  solar  energy  collecting  panels  are  utilized,  and  if  the  center  of 
solar  radiation  pressure  does  not  coincide  with  the  mass  center  of  the  vehicle.  The 
torque  produced  would  be  periodic;  the  period  being  the  time  of  one  orbital  revolution. 

2 

For  a  satellite  employing  flat  solar  energy  collecting  panels,  100  feet 
in  area,  with  a  2.  5  inch  offset  from  the  mass  center,  a  torque  of  approximately 
60  dyne-cm  can  be  produced.1 
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(c) Earth's  Magnetic  Field 

There  are  two  primary  sources  of  torque  resulting  from  the  interaction 
with  the  earth*  s  magnetic  field.  The  first  is  produced  by  interaction  with  the  current 
loops  and  ferromagnetic  material  in  the  vehicle  and  the  second  is  produced  by  eddy 
currents  set  up  by  motion  of  electrically  conducting  parts  of  the  vehicle  relative  to 
the  magnetic  field.  The  latter  becomes  more  prevalent  as  the  angular  velocity  of  the 
vehicle  increases.  The  torque  produced  depends  heavily  on  the  design,  altitude  and 
latitude  of  the  vehicle.  If  large  current  loops  are  avoided  and  nonmagnetic  materials 
are  used  whenever  possible,  the  torque  can  be  kept  to  the  order  of  10  dyne -cm.  As 
in  the  previous  case,  the  torque  produced  is  periodic 

(d)  Internal  Moving  Parts 

There  is  an  angular  momentum  associated  with  every  moving  part  in  the 
vehicle.  The  variation  of  angular  velocity  of  tape  recorder  reels  or  any  rotating 
machinery  causes  a  change  of  angular  momentum  with  respect  to  inertial  space.  This 
in  turn  produces  a  reaction  torque  which  is  exerted  on  the  vehicle. 

(e)  Meteorite  Impacts 

Meteorite  impacts  will  produce  a  torque  on  the  vehicle-  The  magnitude 
of  this  torque  will  vary  in  accordance  with  the  kinetic  energy  of  the  meteorites  and 
the  location  of  impact.  The  time  average  of  the  resultant  torque  is  of  statistical 
nature . 

(f)  Gravity  Gradient  (Differential  Gravity) 

A  portion  of  the  vehicle  a  small  distance  nearer  the  earth  than  another 
part  will  experience  a  slightly  larger  gravitational  force.  Depending  on  the  vehicle 
attitude  and  geometry,  a  resultant  torque  might  be  produced.  The  magnitude  of  this 
torque  will  decrease  as  the  altitude  of  the  satellite  is  increased. 

There  are  a  number  of  basic  methods  which  can  be  utilized  to  obtain  attitude 
control.  These  include: 

(i)  Mass  Expulsion 

This  method  utilizes  familiar  rocket  principles  for  the  generation  of  attitude 
control  torques.  Reaction  jets  are  provided  and  a  propellant  is  stored  within  the 
vehicle.  This  type  of  system  is  good  for  intermittant  large  angle  corrections. 
Applications  of  this  system  are  limited  since  the  propellant  is  exhaustable. 

(ii)  Momentum  Exchange  Within  the  Vehicle 

This  method  utilizes  the  motion  of  internal  parts  to  provide  torques.  Extern¬ 
ally  applied  torques  can  be  countered  by  changing  the  angular  momentum  of  internally 
mounted  gyroscopes  or  inertia  wheels. 
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(iii)  Interaction  with  the  Environment 

This  method  utilizes  prevailing  environmental  conditions  to  provide  torques. 
A  classic  example  is  the  stabilization  of  aircraft  using  atmospheric  drag.  Satellites 
operating  above  an  altitude  of  a  few  hundred  miles  do  not  experience  enough  atmos¬ 
pheric  drag  to  be  utilized  effectively  as  a  means  of  stabilization.  However,  since 
gravity  gradient  torques  tend  to  align  the  axis  of  least  inertia  along  the  geocentric 
vertical,  it  can  effectively  be  utilized  in  controlling  the  orientation  of  the  vehicle. 

The  system  with  which  this  report  is  concerned,  utilizes  gravity  gradient 
torques  to  control  the  vehicled  orientation  and  gyroscopes  to  provide  vehicle 
damping.  Several  different  gyroscope  configurations,  utilizing  single  degree  of 
freedom  gyroscopes  are  considered,  and  for  these  configurations  the  equations 
governing  the  statics  and  dynamics  of  the  system  are  systematically  developed. 

This  report  is  subdivided  in  the  following  manner:  In  Chapter  II  the 
dynamics  of  a  rigid  body  orbiting  the  earth  are  developed  in  terms  of  the  attitude 
deviation  angles  and  the  parameters  of  the  orbit.  Chapter  III  presents  a 
simplified  analysis  of  gyroscope  damping.  In  Chapter  IV,  the  dynamics  of  a  gyro¬ 
scope  constrained  within  a  satellite  are  developed;  the  equations  are  obtained  in 
terms  of  satellite  attitude  deviation  angles,  orbit  parameters,  and  the  gyroscope 
precession  angle.  The  equation  of  constraint,  relating  the  gyroscope  precession 
angle  to  the  satellite  deviation  angles  is  also  developed.  In  Chapter  V,  steady  state 
requirements  for  an  orbiting  vehicle  with  internally  mounted  gyroscopes  are 
considered.  The  equations  obtained  in  Chapters  II,  IV  and  V  are  combined  in 
Chapter  VI  to  obtain  the  dynamics  of  an  orbiting  vehicle  for  various  gyroscope 
configurations . 


See  reference  14  . 
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Chapter  II  -  SATELLITE  ATTITUDE  DYNAMICS 

A)  Reference  Coordinate  System 

In  order  to  develop  a  set  of  equations  suitable  for  the  dynamic  analysis 
of  a  given  system,  it  is  necessary  to  introduce  a  coordinate  reference  for  that  system- 
For  the  problem  of  attitude  control,  it  is  necessary  to  define  a  coordinate  reference 
system  from  which  attitude  deviations  can  be  measured.  For  geocentric  attitude 
control  it  is  desirable  for  the  reference  coordinate  system  to  be  chosen  coincident 
with  the  desired  orientation  of  the  satellite-  The  desired  orientation  of  the  vehicle 
is  determined  by  the  instantaneous  geocentric  vertical,  therefore  the  reference 
coordinate  system  is  not  fixed  and  will  in  general  be  a  variable  angular  velocity 
system.  With  this  choice  of  reference  coordinate  system,it  will  be  possible  to 
linearize  the  resulting  equations  for  small  angular  deviations,  (cos  0=1,  sin  0  =6  )• 
Indeed,  for  large  angular  deviations  these  approximations  cannot  be  made  and  the 
resulting  equations  of  motion  will  contain  nonlinear  trigonometric  terms. 

The  reference  coordinate  system  used  (which  will  be  denoted  by 
X  Y,  Z)  is  shown  for  an  arbitrary  position  in  the  orbit-  (See  Figure  2-1.)  The  Z  axis 
is  taken  along  the  instantaneous  geocentric  vertical,  positive  outward.  The  X  axis 
is  chosen  tangential  to  the  orbit  and  in  the  general  direction  of  the  orbital  velocity  of 
the  vehicle.  The  Y  axis  is  chosen  perpendicular  to  the  plane  of  the  orbit,  in  the 
direction  of  the  orbital  momentum,  such  that  a  right-handed  system  is  formed.  This 
system  is  both  translating  and  rotating  with  respect  to  inertial  space.  It  should  be 
noted  that  the  X  axis  does  not  coincide  with  the  direction  of  the  velocity  vector  of  the 
vehicle  for  an  elliptical  orbit- 

B)  Body  Coordinate  System 

It  is  also  necessary  to  define  a  coordinate  system  which  is  fixed  with 
respect  to  the  satellite  and  which  is  to  be  controlled  to  coincide  with  the  reference 
coordinate  system-  It  is  advantageous  to  choose  the  body  fixed  coordinate  system 
to  coincide  with  the  principal  axes  of  the  vehicle.  This  eliminates  all  product  of 
inertia  terms  and  hence  simplifies  the  differential  equations  of  motion.  (The  inertia 
tensor  is  then  diagonal;  see  Appendix  A).  This  choice  of  body  fixed  axes  reduces  the 
generality  of  the  results, since  it  now  requires  that  the  control  axes  be  principal  axes 
of  the  body.  Since  most  satellite  vehicles  have  some  sort  of  symmetry,  this 
reduction  in  generality  is  not  serious.  It  is  further  assumed  that  the  satellite  is  a 
rigid  body-  This  is  not  strictly  true  due  to  motion  of  its  internal  components  (such 
as  tape  recorders  and  gyroscopes).  If  the  mass  of  such  internal  components  is  small 
in  comparison  with  the  overall  vehicle  mass,  then  these  effects  can  be  neglected- 
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Fig.  2-1  Satellite  Reference  Coordinate  System 
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It  can  therefore  be  assumed  that  the  principal  axes  are  fixed  in  the  satellite.  The 
coordinate  system  formed  by  the  body  principal  axes  is  denoted  by  X^,  Y^,  and  Z^  • 

C)  The  Eulerian  Angles  and  Transformation  Matrix 

The  body  principal  axes  can  be  related  to  the  reference  coordinate  system 

by  an  array  of  Eulerian  angles;  6,0,6  These  angles  are  defined  by  successive 

x  y  z 

rotations  of  the  body  axes  through  the  angles  6^,  6^  and  0^  about  the  X^»  Z^ 
body  axes  respectively  which  rotate  the  vehicle  from  its  reference  system  into  its 
present  position\(see  Fig.  2-2)*  These  angles  will  be  taken  as  positive  when  the 
sense  of  the  rotation  is  equivalent  to  that  of  a  right-handed  screw  advancing  in  the 
direction  of  the  axis  of  rotation.  Figure  2-2  shows  the  three  successive  rotations 
and  their  corresponding  transformation  matrices.  Viewing  the  position  vectors  as 
column  matrices,  the  transformation  matrix  gives  the  new  coordinate  system  as  the 
product  of  the  transformation  matrix  and  the  old  coordinate  system. 


(2.1) 


The  overall  transformation  matrix  relating  the  components  of  a  vector  in  the 
X  Y  Z  system  to  the  componets  of  the  same  vector  in  the  system  is  given  by 


A  = 


cos  6  cos  0 

y  z 


-  cos  6  sin  0 

y  * 


sin  6 


A  =  <v (AJ 

[sin  0  sin  0  cos  6 
y  x  z 

+  cos  0  sin  0  1 
X  zj 

[-  sin  0  sin  0  sin  0 
l  x  y  z 

+  cos  0  cos  0  1 

X  %\ 

-  sin  6  cos  6 

x  y 


(2.2) 


[■ 


cos  6  sin  6  cos  6 
x  y  z 

4-sin  0  sin  0  1 
x  s  J 

cos  6  sin  0  sin  0 
x  j  a 

+•  sin  0  cos  0 
X  z 

cos  6  cos  0 

*  y 


(2.3) 


The  matrix  multiplications  are  not  commutative  and  therefore  must  be  performed  in 
the  order  indicated.  The  physical  significance  of  the  non-commutability  of  the  matrix 
produces  that  the  order  of  the  rotations  cannot  be  interchanged  and  must  be 
performed  as  indicated. 

If  the  deviation  angles  are  small  (sin  0=6,  cos  6  =  1)  the  transformation 
matrix  (2.  3)  bee  omes 

'  \  e_  e  +  0  -0+60 


.0 

z 

®y 


•  e  +  e 

X  y  z 

-  6  6  0  +1 

x  y  z 


-  e* 


(2.4) 
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Upon  further  assuming  that  the  product  of  two  small  angles  can  be  neglected 
compared  to  a  single  angle,  one  obtains 


A  = 


/ 1 

-e 


e 


(2.5) 


-e 


\  y  x 

D)  Development  of  Equations  of  Motion  of  the  Body 

The  equations  of  motion  of  a  rigid  body  in  three  dimensional  space  is 
given  by  Euler’s  equations.  (A24t  A25,  A26)  .  A  derivation  of  these  equations  is 
given  in  appendix  A  . 


T.  =  I  c*)v  +  (I 

iu  X,  b.  z 

b  b  ib 

T.  =  I  +  (I 

Jb  yb 


-  I  )  W 
b  yb  j 


b.  wb. 
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U3B  Wb 
b  k. 


*b 

yb 


'V  \ 


T,_  =  L,  wu  +  (I„  -  I„  )  wb 


Jb' 


(2.6) 

(2.7) 

(2.8) 


D  *'D  D 

where  I  ,  I  I  are  the  moments  of  inertia  about  the  X,  ,  Yu,  Zu  principal  axes 
x,  y,  z,  b  b  b 


of  the  body. 


cj  ,  <*)  ♦  are  the  X,  »  Y,,  Z,  components  of  the  angular  velocity 

T  1  b  b  b 


b^b  kh 

°of  the  body  relative  to  inertial  space. 

T.  ,  T 


Jb1 


are  the  X^»  Y^»  Z^  components  of  the  torque  taken  about 


iv-  W  U  W 

the  mass  center  of  the  body. 


In  order  to  adapt  Euler' s  equations  to  the  problem  of  attitude  control, 
it  is  necessary  to  express  the  angular  velocities  9  wb  9  wb  ^  anc*  tlieir  derivatives 

ib  jb  kb 

)  in  terms  of  the  deviation  angles  (G  ,  0  ,  0  )  and  the  angular  velocity 

D .  D.  D«  X  V  Z 

*b  ^  ^ 

of  the  reference  system.  The  angular  velocity  of  the  body  with  respect  to  inertial  space 
is  given  as  the  vector  sum  of  the  angular  velocity  of  the  vehicle  with  respect  to  the 
reference  frame  and  the  angular  velocity  of  the  reference  frame  relative  to  inertial  space. 


w./inertial  =  w,  //reference  +  0)  reference  /  inertial  (2.9) 

b /  spice  b  /  frame  frame  /  space  1  1 

The  angular  velocity  of  the  body  relative  to  the  reference  frame^is  given  by 

“b /r*r£*”ence  =  ®x  *b  +  ey  jb  +  ez  kb  <2.  10) 


9 

The  angular  velocity  of  the  reference  frame  with  respect  to  inertial  space  is  given  by 


ref  /inertial 
space 


n  i  +  Q  j  +  ft  k 

X  y  J  z 


(2. 11) 


where  are  the  components  of  the  angular  velocity  of  the  reference  frame 

with  respect  to  inertial  space  along  the  X,  Y,  Z  axes  of  the 
reference  system  respectively 


Using  the  transformation  matrix  (2.5),  the  angular  velocity  of  the  reference  system 
can  be  expressed  in  terms  of  the  body  coordinates. 


x\ 

l 

I1  K  -\\ 

!°A 

X 

\xl 

11 

-e  i  e 

z  X 

l  e  -e  i  / 

V y  x  / 

(2.  12) 


OJ 


ref 


+  e  n 

z  y 


e  n 
Y  2 


K 


+  f  -  e  n  +n  +  e  n  1  j. 

[  z  x  y  x  zj  b 

+  [e  n  -  e  n  +  «  1  k. 

[  y  x  x  y  zj  b 


(2.  13) 


Addition  of  Equations  (2.10)  and  (2.13)  ,  and  separation  into  components  yields  the 
following  equations  for  the  components  of  the  angular  velocities  and  their  derivatives 
along  the  body  axes. 


cj  ~  0  -ff2  +  e  n  -Go 

b.  x  x  z  v  v  z 

*b  . 


’b. 


Jb 


=  e  +rc  +  en  -  e  n 

y  y  x  z  z  x 


=e+n+en-en 


y  x  X  y 


=  e+n_  +  en+en  -  e  n 

x  z  y  z  y  y 


=e+n+en4en-en 

y  y  xz  xz  zx 


=  e  +  n  +  en  +e«  -  e  n 

z  z  yx  yx  xy 


e  fi 
y  z 


e  n 

Z  X 


e  n 
*  y 


(2.  14) 
(2.15) 
(2.  16) 
(2.  17) 

(2.  18) 
(2.19) 
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Substituting  these  results  into  Euler's  equations  (2.6,  2.7,  2.8)  yields: 

t.  =  i  [e  +  h  +  eii  +  e  n  -  e  n  -  e  n  1 

l,  x,  1  x  x  zy  zy  yz  yzJ 


fee  +eo  +een  +een  +en 

L  y  z  yz  xyy  y  y  x  zy 


+  (I  -  I  ) 
zb  ^b 

+  n  n  -en2+enn  -  e  e  n  -  e  n  n 

yz  x  y  y  x  y  zzx  zxz 

+  eenn  -  e  e  o2  +  e  e  n  +  e  n  2  -  e  2  n  « 

xzxy  yzx  xzz  xz  xyz 

+  e  e  q  n  1 

X  y  x  zj 

t  =i  fe+n-en-en+en+6n|| 
jk  yK  u  y  y  z  x  z  x  x  z  x  zJi 


Jb  ^b 

+  (I 


-  i 

b  zb 


)fee  +  eo  -een  h  en  +  e  n 

[  x  z  xz  xxy  y  x  x  zx 


+  n  n  -enn  +en+een  +  e  n  n 

xz  xxy  y  x  zzy  zyz 

-een2+eenn  -  een  -  en2-e2nn 

xzy  yzxy  yzz  y  ®  yxz 

+  o  e 
*  y 


n  n  ) 

y  zjl 


t  =i  [e-n-en-en+en  +  e in  1 

k.  zk*-z  z  xy  X  y  yx  yxj 


)  fe’  e  +  e 

/ 1  x  y 


en  -een  +een  +en 


+  (i  -  i 
yb  J 

+  n  n  -en2+enn  +een  +en2 

xy  zx  xxz  zyy  zy 

-e2nn  +  eenn  -  een  -  enn 

z  xy  xzyz  yyz  yyz 

+  eenn  -een2! 

.  — r  z  j 


y  z  x  z 


x  y 


(2.20) 


(2.21) 


(2.22) 


where  T.  »  T.  •  T.  are  the  components  of  torque  about  the  center  of  gravity  of  the 
lb  Jb  o 


vehicle  alon^  X^#  and  respectively. 


The  above  equations  are  valid  under  the  following  assumptions: 

a)  Inertia  of  internal  moving  parts  of  the  body  are  negligible.  - 
This  assumption  is  necessary  to  allow  the  system  to  be  treated  as  a  rigid  body. 

b)  The  fixed  body  axes  coincide  with  the  principal  axes.  - 
This  assumption  made  the  enertia  tensor  diagonal  (see  eq.  A.  21). 

c)  Small  angular  deviations.  - 

This  assumption  is  utilized  in  the  de rivation  of  the  transformation  matrix  (2.  5)  . 

In  order  to  specialize  these  equations  for  the  problem  of  attitude  control^  it 
is  necessary  to  obtain  an  expression  for  the  components  of  the  reference  system  angular 
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velocity  (ftx#  ft^,  ft^).  If  the  earth  were  a  perfect  homogeneous  sphere^*  and  if  the 
only  force  acting  on  the  satellite  was  the  earth’s  gravity,  the  satellite  orbit  would  be 
an  ellipse  whose  spatial  plane  and  major  axis  would  be  fixed  in  inertial  space.  Under 
these  conditions,  ftx  and  f 2^  would  be  zero  and  ft  would  be  equal  to  the  orbital 
angular  velocity  (sometimes  referred  to  as  the  pitch  over  rate)  predicted  by  Kepler' s 
laws  of  planetary  motion.  There  are  several  perturbations  which  act  on  the  orbit  and 
change  ftx»  ft^  and  ft^  from  these  idealized  values .  They  are  listed  below: 

a)  The  earth  is  not  a  perfect  sphere  but  is  closer  to  an  oblate  spheroid. 

b)  The  earth  is  not  homogeneous. 

c)  Atmospheric  drag. 

d)  Gravitational  attraction  of  the  moon,  sun  and  other  planets. 

e)  Electromagnetic  forces  caused  by  the  satellite's  motion  through 
the  earth's  magnetic  field. 


f)  Collision  with  meteorites  and  other  particles. 


The  changes  in  angular  velocity  produced  by  the  perturbations  discussed  above 
are  several  orders  of  magnitude  less  than  the  orbital  angular  velocity  which  is  in  the 
order  of  10  rad/  sec.  These  perturbations  therefore  have  a  negligible  effect  in  the 
equations  of  motion  when  compared  to  terms  containing  ft  .  As  a  result,  ft  and 
a rC  cons^erec*  to  be  zero  in  the  equations  of  motion.  It  is  further  assumed  that 
there  is  no  time  average  force  acting  on  the  vehicley since  such  a  force  would  cause  the 
plane  of  the  orbit  to  become  time  varying. 

Setting  ftx  =  ftx  =  ftz  a  ft^  =  0  in  the  torque  equations  (2.20,  2.21,  2.22) 

yield : 


T. 

lb 


+  e  ft 

z  y 


+  0  ft 
z  y 


n  -  e  «  2 
y  * 


e  +  e  e  n 

z  x  y  y 


] 


(*•23) 


i  [  e  +  n  1  +  (i  -  i  )  fe  e  -  e  e‘  n 

Vb  l  y  y  J  xb  *b  lxi  x  x  y 


+  e  e  n  -  e  e  n  2 1 

z  z  y  x  z  y  J 


i„  f  -  ®  n  *  6  n  1  +  (i  -  i  )  f  e  e  +  e*  n 

V  z  X  y  X  yj  yb  XK'  L  X  y  x  • 


+  o  6  «  +e  n  21 

z  y  y  z  y  J 


(2.24) 


(2.25) 
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Further  simplification  of  Eqs.  (2.23),  (2.  24),  (2.  25)  can  be  obtained  if  the 
order  of  magnitude  of  the  various  terms  is  estimated. 

1)  fly  is  taken  to  be  10  ^  rad/sec.  This  corresponds  roughly  to  a  ninety 
minute,  400  mile  high  circular  orbit. 

2)  An  estimate  of  the  angular  acceleration  can  be  obtained  from  a  simple  example. 
Consider  a  1000  pound  drum  shaped  satellite  (see  Fig.  2-3)  which  is  4  feet  in  length 
and  4  feet  in  diameter 


Fig.  2-3  Drum  Shaped  Satellite 


For  a  cylindrical  satellite 


rx  =  I  =  +3R2) 

xb  yb 

I  -imR2 

zb  2 


(2. 26) 
(2.27) 


Where  I  ,  I  and  I  are  the  polar  moments  of  inertia  about  the  X.  Y  and  Z 
b  "b  zb  b  nt>  b 

axes  of  the  body. 


m  is  the  mass  of  the  satellite  in  slugs. 
I  is  the  length  of  the  satellite  in  feet. 

R  is  the  radius  of  the  satellite  in  feet. 


Calculation  of  the  moments  of  inertia  yields: 
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I  «1  -72.5  slug  -  ft.2  (2-  28) 

*b  Vb 

I  =62.2  slug  -  ft.2  .  (2  29) 

ZV 


It  will  be  assumed  that  a  sinusoidal  disturbance  torque  of  amplitude 
80  dyne  -  cm  is  applied  to  the  vehicle  with  a  radian  frequency  of  O  (10~3  rad/ sec) 
about  its  axis  . 


80  sin  10“3t  dyne -cm 


(2.30) 


Converting  to  ft.  -  lbs.  yields: 

Td  =  5.9  x  10*6  sin  10*3t  ft  -  lbs. 
*b 


(2.31) 


Assume  for  the  purpose  of  estimation  that  the  elementary  equation 

applies. 

Therefore  the  angular  acceleration  is  given  by 


0.94  x  10'7  sin  10'3t  — 
sec 


(2.32) 


(2.33) 


The  maximum  value  of  the  angular  acceleration  about  all  the  axes  is  seen  to  be  about 
IQ'7  rad/  sec^ 
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-7  rad 

2 

sec 


(2- 34) 


3)  An  estimate  of  the  angular  velocities  and  angular  deviations  can  be  obtained 

by  integration.  Integrating  to  obtain  the  angular  velocity  and  deviation,  one  obtains 


&  =  -  10'4  cos  10"3t  (2. 35) 

6  =  -  10*  1  sin  10‘3t  (2.  36) 


Therefore  the  estimate  for  the  maximum  value  of  the  angular  velocity  will  be 
10  rad/ sec. 

0=6=0=  10~*  rad/  sec 
x  y  z 


(2.37) 


and  the  estimate  for  the  maximum  value  of  the  angular  deviation  will  be  0.  1  rad 


0  =  G  =  G  =  0.  1  rad 

x  y  z 


(2.  38) 
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The  effect  of  damping  on  the  vehicle  will  cause  these  values  to  be  smaller  than  those 
used  for  the  estimations. 


Substituting  the  estimated  values  into  equations  (2.23),  (2.24)  and  (2.25) 
and  deleting  the  negligible  terms,  (a  term  will  be  considered  negligible  if  it  is  at 
least  one  order  of  magnitude  smaller  than  the  other  terms),  yields  the  final  set  of 
equations  of  motion  for  a  body  in  a  general  elliptical  and  planar  orbit. 


T  21  (2'39) 

t.  =  i  fd'+efl+enl  +  d  -i)6n-en 
•b  *bl*  *  *  r  vl  S  *  v  J  (2  40) 

\  =  Vb  [  \  *  riy]  *  \  '  V  [#.  9.  -  V*“y  * VV  -  6x»A2] 


T,  =  i  fe-eft-enl  +  (i  -i  )  [e  n  +  e  n  2*[  (2.4i) 

\  zb  I  z  x  y  x  yJ  yb  xb  L x  y  z  y  J 

E)  Inclusion  of  the  Gravity  Gradient  Torque 

It  is  essential  to  understand  that  the  left  side  of  equations  (2.  39),  (2.40)  and 
(2.41)  represent  the  components  of  the  total  external  torques  applied  to  the  body. 


Fig.  2-4  Mechanism  of  Differential  Gravity  Torque 


IS 


In  general,  the  total  external  torque  applied  to  the  vehicle  is  composed  of  disturbance 
torques  and  control  torques.  It  is  advantageous  to  explicitly  include  the  expressions 
for  the  control  torques  so  that  the  only  torque  appearing  on  the  left  side  of  the 
equations  are  the  true  disturbance  torques.  , 

Differential  gravity  is  a  possible  source  of  a  control  torque-  The  mechanism, 
by  which  differential  gravity  can  lie  utilized  to  provide  a  control  torque,  can  be 
qualitatively  demonstrated  with  the  aid  of  the  following  simple  example: 

Consider  the  dumbell  shaped  vehicle  as  shown  in  Fig.  2-4.  The  difference  in 
force  acting  at  itsendswill  cause  a  small  torque  tending  to  align  the  longitudinal 
axis  of  the  vehicle  along  the  geocentric  vertical.  As  the  satellite  orbits  the  earth  a 
torque  is  thereby  produced  tending  to  keep  the  satellite  so  aligned.  The  positions  of 
stable  and  un- stable  equilibrium  are  also  shown  in  Fig.  2-4.  It  is  of  cardinal 
importance  to  realize  that  differential  gravity  alone  can  be  used  as  a  mechanism  for 
keeping  an  axis  of  the  satellite  aligned  in  the  general  direction  of  the  earth.  Without 
damping  however,  the  poles  of  the  system  will  lie  on  the  imaginary  axis  and  therefore 
a  disturbance  torque  will  cause  the  vehicle  to  oscillate* 

A  derivation  for  the  differential  gravity  torque,  for  a  sufficiently  general 
body,  is  given  in  Appendix  C.  The  results  of  the  derivation  are  repeated  here  for 
c  onvenience.  For  an  elliptical  orbit 


3  g  R 

V^rr- 


i,  e  (i  -  i  )  +  j.  e  (i  -  i  ) 

b  x  zb  V  b  y  zt  xb 


(C22) 


and  for  a  circular  orbit 


T  ,  = 

d  y 

g  y 


i  0  (I  -  I  )  +  j,  0  (I  -  I  ) 

b  x  zb  Vb  b  V  zb  xb 


(C23) 


It  is  immediately  seen  from  these  equations  that  the  component  of  the  gravity  gradient 
torque  about  the  axis  of  the  body  is  zero.  Therefore,  differential  gravity  cannot 
l)e  used  to  prevent  drift  about  the  Z^body  axis.  It  is  seen  from  either  equation  (C22) 
or  (C2  3)  that  differential  gravity  introduces  a  component  of  torque  on  both  the  X^and 
axes  of  the  body  that  is  proportional  to  the  displacement  angle  about  that  axis. 

Since  the  behavior  of  the  differential  gravity  torque  on  both  the  and  Y^  axes  are 
similar,  only  the  component  along  the  Y^  axis  will  be  considered. 
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Two  situations  are  of  interest: 


1)  I  <  I 
xb  zb 


2)  I  <  I 


zb  xb 


Consider  case  (1) 


+  3  0  II 


(2.42) 

(2.43) 


(2.44) 


It  is  seen  from  equation  (2.44)  that  the  differential  gravity  torque  is  in  the  same 
direction  as  the  angular  displacement.  This  clearly  represents  an  unstable 
configuration  since  once  a  disturbance  torque  produces  a  small  angular  displacement, 
differential  gravity  will  cause  this  angle  to  increase. 


Now  consider  case  (2) 


(2.45) 


In  this  case  differential  gravity  is  in  the  opposite  direction  of  the  angular  displace¬ 
ment.  Hence  if  a  disturbance  torque  produces  an  angular  displacement,  differential 
gravity  will  tend  to  restore  the  vehicle  to  its  unperturbed  position.  For  this  configu¬ 
ration,  differential  gravity  introduces  position  feedback  much  in  the  same  manner 
as  a  body  coupled  to  a  spring.  Therefore,  it  is  of  utmost  importance  that  the  moment 
of  inertia  about  the  axis  of  the  body  be  less  than  the  moments  of  inertia  about 

both  the  X,  and  Y,  axes  of  the  body, 
b  b 


From  equation  (C  22)  it  is  seen  that  the  total  external  torque  on  the  body  is 


given  by 


3  g  R 

t.  »  t,  +  — V-e  e  (i  - 1  ) 

*b  di.  P  3  x  zb  yb 
b  r  o 


3  g  R 

Ti  *  Td  +— V*  6y(Iz  -!xj 

•>b  jb  poJ  b  b 


(2.46) 


(2.47) 


(2'48> 

where  T^  is  .the  total  torque  on  the  body  less  the  differential  gravity  torque. 
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Substituting  equations  (2.46),  (2.47)  and  (2.48)  into  equations  (2.  39),  (2-40)  and 
(2.4l),  and  bringing  the  differential  gravity  terms  over  to  the  right  side  of  the 
'■'luations  yields: 


r  =  i  +  0  fi  +  G  ft]  +(I  -I  )  [6  fl 

d  xb  l  x  z  y  z  yj  zfe  yb  L  z  V 


2  3gs  Re  1 
-  0  (U  6  + — - — 2  )| 

x  y  «  3  J 


(2.49) 


e  e  n  +  e  e  n 

x  x  y  z  z  y 


t  .  =  i  +  n  I  +  (i  -  i  )  [6  6  - 

dj.  yb  l  y  yJ  ^  v  Lx  2 

-  e  e  n  2  +_! — *  e  l. 

x  z  y  p  3  y  J 

Ko 

t.  =  i  [8  -  e  n  -  e  o  1  +  (i  - 1  )  [fi  n  +  e  n  21  (2  s 

\  zb  l  z  X  y  X  yj  yfe  xb'  [ x  y  z  y  \ 

The  corresponding  equations  of  motion  for  a  circular  orbit  can  be  obtained 


(2. 50) 


1) 


by  setting  =  0  and 


3gs  Re  2 

— 2 — i  =  3  n  c 


t  .  »  i  f8  +  e  n  1  +  (i  -  i  )  [6  n  -  4  e  n  2 1 

d.  xb  l  x  z  yJ  Z  yb  L  z  y  x  y  ] 
b  b 


t  ,  ~  i  8  +  (i  -i  )  r§ 


6-eonteen 

xz  xxy  zzy 


0  0  fl  +36 

x  z  y  y 


vi 


t  .  =  i  fe  -  e  n  1  +  (i  - 1  )  n  +  e  n  21 

\  zb 1 2  x  yJ  yb  xb  L x  y  2  y  i 


(C.16) 

(2.52) 

(2.53) 

(2.54) 


Equations  (2.  50)  and  (2.  53)  can  again  be  simplified  if  order  of  magnitude  approximations 
are  considered.  Using  the  estimated  values  given  in  equations  (2.  37)  and  (2.38)  it 
is  seen  that  equation  (2.50)  simplifies  to  2 

(2.55) 


6  g  R 

t.  =i  fe  +61  +3(1  .  i  )  JLe  e 

dj  yb  l  y  yJ  x*  zk  „ 


V,  6, 

b  bp 


and  equation  (2.53)  simplies  to 


T ,  =  I  6  +  3(1  -  I  )  6  0 

djb  yb  y  xb  zb  y  y 


2 


(2.  56) 
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It  should  be  noted  from  equations  (2*  49),  (2.  51)  and  (2.  55)  that  the 
expressions  for  and  are  linear,  cross  coupled,  time  varying  differential 


lb 


equations  that  depend  on  0^  and 


derivatives.  Therefore,  0 


6  and  their  derivatives  but  not  on  G  and  its 
z  y 

and  0  can  be  determined  from  the  disturbance  torques 


about  the  X^  and  Z^  axis  without  involving  the  disturbance  torque  about  the  axis. 

Equation  (2.55)  is  a  linear,  time  varying  differential  equation  that  only  involves  0^ 
and  its  derivatives  and  does  not  depend  upon  the  angular  deviations  about  the  other  two 


As  an  aid  to  obtaining  an  insight  into  equations  (2.49),  (2.51)  and  (2.  55),  they 
are  put  into  the  following  form,  and  from  this  form,  a  signal  flow  graph  can  be 
easily  constructed,  (see  Fig.  2-5  and  Fig.  2-6) 

T, 


d. 

i, 


e  = _ k  .  e  h 


. . .  .(v_y 


z  y  z  y 


,  2  38s  Re 

e  n  -e  n  c  + — - — 
z  y  *  y  „  3 


(2.57) 


e_  Jl?. n  -  3^) 


y  |  y 

*b 


ly 

Yb  '  0 


0  = 
z 


+  G 


»  *  (yb  *v)  r. 

n  +e  n  - 2/1  e  n  +  e  n 

y  *  y  :  y  *  y 


(2.58) 


(2.59) 


and 


Figure  2-5  relates  the  output  angles  0^  and  G ^  to  the  torques  about  the 
axes,  and  Fig.  2-6  relates  6^  to  the  torque  about  the  Y^  axis. 


With  reference  to  Fig.  2-5,  paths  P^  and  represent  position  feedback 

terms  due  to  the  orbital  angular  velocity  of  the  vehicle.  If  I  >  I  and  I  >  I  . 

*b  yb  V 

the  feedback  produced  by  these  paths  will  be  such  as  to  oppose  a  displacement 
produced  by  a  disturbance  torque.  Paths  P^  and  P^  are  also  due  to  the  orbital 
angular  velocity  of  the  satellite  about  the  earth.  These  paths  introduce  feedback 
depending  on  the  error  rates  and  G^)  .  Unfortunately  these  terms  do  not 

introduce  damping  because  they  appear  in  the  wrong  equations.  For  an  elliptical 
orbit,  the  angular  momentum  of  the  vehicle  becomes  time  varying.  The  derivative 

9 

of  the  orbital  rate  (0^)  gives  rise  to  c^oss  coupling  between,  the  and  Z^  axes 
as  shown  by  paths  and  P^.  With  reference  to  both  Figs.  2-5  and  2-6,  paths  P^ 
and  Pg  are  seen  to  introduce  position  feedback  on  the  and  Y^  axes.  These  are 
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I 


Elliptical  Orbit 
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the  result  of  differential  gravity  and,  as  was  previously  shown,  the  feedback  will  be 
negative  about  both  axes  if  the  moment  of  inertia  of  the  body  about  its  Z^axis  is 
smaller  than  the  moments  of  inertia  about  the  other  two  axes.  The  corresponding 
flow  graphs  for  a  circular  orbit  can  be  obtained  from  Figs.  2-5  and  2-6  by  netting 
fly  equal  to  zero  and 


(C.16) 


It  is  seen  from  the  flow  graphs  that  position  feedback  terms,  for  a  Suitable 

choice  of  the  moments  of  inertia,  (I  >  I  >  I  )  can  act  as  spring  type  restoring 

yb  *b  zb 

torques  about  each  of  the  axes  of  the  body.  The  effective  spring  constant,  defined 
as  the  ratio  of  restoring  torque  to  angular  displacement,  is  given  below  for  each  of 
the  axes. 


Kx  =v  -  h  >. 
7  b  b 


2  3gs  Re 

<  +  -rr- 


(2.60) 


K  =  3  (I  -  I  )  .8s.  ,f.e 
y  Xb  zb  p  3 


(2.61) 


K  (I  -  I  )  12 
Z  yb  Xb  y 


(2.62) 


It  is  of  interest  to  obtain  some  appreciation,  as  to  the  order  of  magnitude  of  these 
terms.  Consider  again  the  drum  satellite  of  Fig.2—3.  Assuming  a  circular  orbit 
and  using  the  moments  of  inertia  calculated  for  this  satellite,  yields. 


K  =  4.11  x  10  ‘  5 
x 

Ky  =  |kx  =  3.11  x  10 

K  =  0 
z 


ft.  -lbs . 
rad 

ft.  -lbs . 
rad 


(2.63) 


(2.64) 
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F)  Solution  of  the  Equations  of  Motion  for  a  Circular  Orbit 

Introducing  Laplace  transform  notation  into  equations  (2.  52),  (2.  54)  and 
(2 . 56)  yields : 


(8)  =  !xj  8"  ex(8)  +  8  ez  (8)  °y  +  (I-  -  )  Is  M8>  «, 


d‘b 


z  y  zb  Vl  z  y 


-46  (s)  n  2 
x  y 


(2.66) 


t.  (s)  =  i  s2  e  (b)  +  3(1  -i  )  e  (s)  n 

djb  yb  y  xb  V  y  y 


(2.67) 


=  I 

s2  0  (s)  -  s  G  (s) 

+  (i 

- 1  ) 

zb 

z  x  y 

yb 

xb  L 

s  0  (s)S7 
x  y 


+  6  (s) 

z  y 


(2.68) 


Solving  equations  (2.66),  (2.67)  and  (2.68)  one  obtains: 


Td.  <8> 


s2  i  +  (i  -  i  )  n  2 

zb  yb  xb  y 


-  X” 


it!  (I 

y  : 


+  I  -  I  ) 

xb  zb  yb 


ey  (s)  =  _ 

s2  n  2  (i  +  i  -i  )2+[b2i  +4(i  -  i  )n  2^ 


(2.69) 


y  xb  zb  yb  I-  ^  yb  zb  y 


s2i  +(i  -i  )n2 

zb  yb  xb  y 


Td.  (8> 


6  (b)  = 

ir  '  ' 


s2 1  +3(1  -i  )n 


(2.70) 


'b  Xb 


zb  y 


6z  (B) 


.  V’ 


82  i  +4(i  -  i  )  n  2 

’S,  yb  zb  y  . 


+  T 


.  (s)[  b  n  (i  + 1  -  i  ) 

V  t  y  xb  zb  ybJ 


2  „  2 


B4,  n  ‘  (I  +  I  -  I  )2  +  s2  I  +  4(1  -  I  )£}  2  b2  I  +  (I  -  I  )0 
y  xb  zb  V  L  xb  yb  zb  yJl  zb  yb  *b  y 


(2.71) 

2l 


i 


In  order  to  better  understand  the  nature  of  the  response  of  the  body  to 
disturbance  torques,  several  examples  will  be  considered. 
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Example  No.  1 

The  moments  of  inertia  about  the  and  axes  are  equal  and  greater  than 
ti  e  moment  of  inertia  about  the  axis.  An  impulsive  torque  of  strength  6  is 
applied  about  the  axis  and  no  torques  are  applied  about  the  other  two  axes. 
Impulsive  torques  are  characteristic  of  meteorite  impacts. 


When  appropriate  substitutions  are  made  in  equation  (2.69)  one  obtains  for 


0  (s) 
x 


e  (s)  = 
x  ' 


(4  1  -31) 

2  ,  n  2  "b  % 

s  +12  - 

y  i 


(2.72) 


Since  I  >  I  the  poles  of  0  (s)  lie  on  the  imaginary  axis  and  the  system  will 
xb  zb  x 


oscillate  with  a  radian  frequency  of 


4  1  -  3  I 


xb  zb 


(2.73) 


If  I  <  I  there  will  be  a  pole  of  0  (s)  in  the  right  half  plane  and  the  angular 
xb  zb  x 

deviation  will  increase  with  time.  For  I  >  I 


0  (t)  =  .  °  sin  «  t  U(t) 

X  I  03  O 

*b  ° 


Clearly  from  equation  (2.70) 
6  (t)  =  0 

y  ' 

From  eqviation  (2.71) 


ez(.)  = 


n  6 
JL 


8fs2I  +n2(4l  -31  )1  s  12  (41  -31 

^  y  '  *b  zb  y  xb  2 

to 


(2  (4  1  -31) 

y  %  % 


s  I  6 

x _ 

7!  +  a  2  (4  I  -31  ) 

*b  y  *b  % 


(2.74) 


(2.75) 


(2.76) 
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Therefore  for  I  >  I 

xb  zb 

6  (t)  =  - - -  [l  -  cos  u  t  ]u(t)  (2.77) 

n  (4  I  -  3  I  )  l  °  J 

Y  xb  zb 

F  rom  equation  (2. 74)  it  is  seen  that  the  angular  deviation  about  the  X,  axis  is 

b 

oscillatory  with  zero  mean,  whereas  from  equation  (2.77)  it  is  seen  that  the  deviation 
about  the  axis  is  also  oscillatory  but  with  a  non-zero  mean.  The  results  of  this 
example  are  summarized  in  Fig.  2-7. 

Example  No.  2 

The  same  satellite  configuration  as  in  Example  1  but  now  the  impulsive  torque 
of  strength  6  is  applied  about  the  axis  and  no  torques  are  applied  about  the 
othe r  two  axes . 


Substitution  into  equations  (2. 69)  and  (2.71)^it  is  clearly  seen  that 


0x(t)  =  Gz(t)  =  0 


(2.78) 


Substitution  into  equation  (2.  70)  yields  for  ©^(s) 


0  (s) 

y 


^b 


,  (IXK  -  V 

82  +  3  — - 5-  n  2 


(2.79) 


Therefore 


0V( t)  a  - sin  cj  jt  U(t) 


y  *  CO  ,1 

lyb 


(2.80) 


where 


w  ,  -  n 

1  y. 


3dx  -  it  ) 

xb  zb 


(2.81) 


yb 


(observe  that  w  <  <*>  ). 

'  1  o 

It  is  seen  from  equation  (2.  80)  that  the  angular  deviation  about  the  Y^  axis 
is  sinusoidal  with  zero  mean.  The  result  of  this  example  is  summarized  in  Fig.  2*8  . 


ANGULAR  DEVIATION 


Fig.  2-8  Response  to  an  Impulsive  Torque  of  Strength  8  Applied  About  the  Axis 
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Example  No.  3 

The  same  satellite  configuration  as  in  Example  1,  but  now  an  impulsive 
torque  of  strength  6  in  Applied  about  the  Z^  axis  and  no  torques  are  applied  about 
the  other  two  axes. 


Substitution  into  equation  (2.69)  yields  for  C^(s) 


e  (s) 


-  ft  6 

_ 3L 


s  s  I  +  ft  (4  I  -  3  I  )]  s  ft  (4  I  -31) 

Xl  V  X  7  '  v  X 

h  ’  K  K  i  *  V 


s  I  6 
x 


vb  'Jb 

T 


ft  (4  I  -31  )  !"s  i  +  ft  2  (4  I  -31  )] 

y  x^  z^  '  x^  y  x^  z^  j 


J 


Therefore 


(2.82) 


f.  (t)  =  - 2 - 

x  S2  (4  I  -31.) 
y  xb  zb 


|cos  u  t  -  1  U(t) 


(2.83) 


Clearly  from  equation  (2.70) 

0  (t)  =  0 

y 

Substitution  into  equation  (2.71)  yields  for  O^s) 


o7(S)  = 


1 

6 

s2I  +4(1  -  I  )  Q  2 

X,  X,  7s\  V 

L  b  b  b  y  J 

s2  I  i 

Zb  1 

s2  I  +  (4  I  -  3  I  )  ft 

Xb  xb  zb  y  ! 

(41  -  31  \ 

j  s2+\A_^)n2 


(2.84) 


I  Y  (2.85) 


4(1  -1)6  4(1  -1)6 

_  b  b  _  _ ^b  b 

s2  I  (41  -  31  )  I  (41  Til  )T 

.  7  V  -7  V  '  7  ' 


zb  xb  zb  zb  xb  %  I  s2 


I41  '3I  \ 

i  Xb  Zb) 

+ - j - 


Taking  the  inverse  Laplace  transform,  one  obtains  l.  b 

r  (4  1  -  I  ^ 

5in  „  t  A  lb  _!w  t 

°  i 


oz(t)  • 


(4  I  -31  )  u 

x,  z,  o 

b  b 


U(t)  (2.86) 


It  is  seen  from  equation  (2.86)  that  when  an  impulsive  torque  is  applied  about  the 

axis  drift  occurs.  Differential  gravity  is  not  available  to  help  stabilize  the  vehicle 

about  that  axis.  Furthermore,  since  it  was  assumed  that  1  =  I  ,  the  gyroscopic 

xb  yb 
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action  due  to  the  orbital  rate  of  the  vehicle  was  not  able  to  produce  position  feedback. 

(see  Fig.  2-5  and  equation  (2.62)).  The  results  of  this  example  are  summarized  in 

Fig.  2-9  .  If  the  moments  of  inertia  of  the  vehicle  were  chosen  such  that 

I  >  I  >  I  ,  drift  would  not  occur  about  any  of  the  vehiclefs  axes  as  a  result  of  an 
yb  xb  zb 


impulsive  disturbance.  Consider  the  following  example. 

inertia  of  the  vehicle  such  that  L,  >  I  >  I  and 

~b  xb  zb 


Choose  the  moments  of 


I  +  I  =  I 
xb  zb  yb 


(2.87) 


2  1 

(For  example:  1^  =  ^  I  ,  1^  =  ^  I  )  .  Equations  (2.  69),  (2.70)  and  (2.71)  degenerate 


to 


ex(8)  =  r2,' 


s2  i  +  (i  -  i  )  n  2 

zb  yb  xb  y 


8  I 


W2.88) 


+  4(1  .1  )n  c  Bc  i  +  ci  -  i  )  n  | 

)  yb  zb  y  JL  zb  yb  *b  y  I 


e  (s) 
y 


“F 

s 


Td  (s) 
Jb 


I 

yb 


+  3(1 


(2.89) 


T  j  (s)|  s2  I  +  4  (I  -  I  )  n  2 


ez{8)  "  777 


’I"’ 


yb  zb  y 


(2.90) 


's  i  +4(1  -I  )Q  ][s  I  +  (I  -I  )  Q 

Xb  Vb  Zb  y  J[  Zb  yb  xb  V 


It  is  apparent  from  these  equations,  that  when  the  moments  of  inertia  of  the  body 
satisfy  equation  (2.87),  the  three  axes  of  the  body  are  completely  uncoupled,  i  e.  a 
disturbance  applied  about  one  axis  will  only  produce  a  displacement  about  that  axis 
and  none  about  the  other  axes.  Also,  it  is  seen  that  no  drift  will  occur  about  any  of 
the  axes  due  to  an  impulsive  disturbance.  The  free  response  consists  of  the  sum  of 
two  sinusoids  of  different  frequency.  This  is  certainly  an  improvement  over  the  drum¬ 
shaped  satellite  discussed  in  Example  1,  2,  and  3. 


We  have  seen  in  this  chapter  that,  with  proper  vehicle  design,  the  impulsive 
response  of  the  system  consists  of  only  oscillatory  terms.  Since  the  system  is  linear, 
if  an  arbitrary  torque  is  applied,  the  output  will  consist  of  terms  of  the  same  form  as 
the  input  plus  the  terms  of  the  impulsive  response.  Hence  for  impulse -like  torques, 
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it  the  Z,  Axis 
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low  level  step  torques  and  sinusoidal  torques,  drift  would  never  occur.  If  damping 
devices  are  incorporated  into  the  system,  (this  can  be  done  passively)  it  should  be 
possible  to  remove  the  oscillatory  components  and  keep  the  vehicle  reasonably  aligned 
with  the  geocentric  vertical-.  It  is  possible  to  utilize  gyroscopes  to  provide  damping. 
In  order  to  better  understand  how  this  can  be  effected,  a  simple  analysis  of 
gyroscope  damping  will  be  presented  in  Chapter  III. 
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Chapter  III  -  A  SIMPLE  ANALYSIS  OF  GYROSCOPE  DAMPING 

A)  The  Mechanism  of  Gyroscope  Damping 

The  mechanism  by  which  gyroscopic  action  can  be  utilized  to  provide 
vehicle  damping  can  be  qualitatively  explained  as  follows*  When  a  disturbance  torque 
is  applied  to  the  vehicle,  angular  motion  of  the  vehicle  with  respect  to  inertial  space 
will  occur.  This  motion  will  cause  relative  motion  between  the  gimbal  of  the  gyroscope 
and  the  body.  If  viscous  damping  at  the  gyroscope  gimbal  is  provided  (this  can  be 
provided  through  electrical  means)  a  damping  term  will  couple  into  the  equations  of 
motion  of  the  body. 


Vehicle  damping  can  either  be  of  direct  or  indirect  nature.  By  indirect 
damping,  it  is  meant  that  vehicle  rotation  will  produce  precession  of  the  gyroscope 
which  causes  vehicle  damping  due  to  relative  motion  between  the  gyroscope  and  the 
vehicle.  If  the  vehicle  rotation  is  about  the  precession  axis  of  the  gyroscope,  direct 
damping  will  occur.  In  this  case  the  relative  motion  between  the  vehicle  and  the 
gyroscope  gimbal  is  brought  about  without  the  mechanism  of  gyroscope  precession. 


B)  Analysis  of  a  Simple  Configuration  Employing  Indirect  Gyroscope  Damping 
Indirect  damping  can  be  demonstrated  with  the  following  example. 
Consider  the  geometry  shown  in  Figure  3-1  .  In  the  geometry  pictured, the  body  is 
constrained  to  rotate  with  respect  to  inertial  space  only  about  its  X  axis  and  the 
gyroscope  is  constrained  to  rotate  relative  to  the  body  about  the  Z  axis  of  the  body. 
Consider  a  disturbance  torque  applied  to  the  body  about  its  X  axis.  This  torque  is 
countered  by  the  inertia  torque  of  the  body,  the  spring  torque^and  the  gyroscopic  torque. 


t  ,  =  i  e  +  k  e  -  e  h  cos  e 

d  x  x  x  g  r  g 


(3.  1) 


Where  1^  is  the  sum  of  the  moments  of  inertia  of  the  body  and  the  gyroscope 
about  the  X  axis . 


K  is  the  spring  constant 

Hr  is  the  angular  momentum  of  the  gyroscope  rotor. 
The  torque  on  the  gyroscope  about  the  Z  axis  is  given  by: 


g 


Where  D  is  the  damping  coefficient  between  the  body  and  the  gyroscope  gimbal. 
This  torque,  (3.  2),  is  countered  by  the  inertia  torque  of  the  gyroscope  and  the 
gyroscopic  torque. 
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Fig.  3-1  Configuration  to  Illustrate  Indirect  Gyroscope  Damping 
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=  -  D0  =  I  G  +  G  H  cos  G 
g  g  g  x  r  g 


(3.3) 


where  1^  is  the  moment  of  inertia  of  the  gyroscope  about  the  Z  axis  of  the  body. 
Equations  (3.  1)^(3.  3)  can  be  linearized  for  small  angular  deviations. 


T,  =  I  G  +KG  -  GH 
d  x  x  x  g  r 


-  D  0  =  I  G+GH 
g  g  g  x  r 


Introducing  Laplace  transform  notation  to  these  linearized  equations  yields: 

Tj(s>  =  i  s2  e  (s)  +  k  e  (s)  -  h  s  e  (s) 

d'  x  x'  9  x'  '  r  g  ' 

-  Ds  6  (s)  =  Is^  G  (s)  +  H  s  G  (s) 
g  g  g  r  x' 

Solving  for  0  (s)  in  equation  (3.7)  in  terms  of  G  (s)  yields; 
x  g 

-  H 

e  (8)  =  — S-  e  (s) 

g'  7  D  +  si  x'  ' 

6  g 

Substituting  into  Equation  (3.6)  yields: 

Td(s)  [d  +  8  Ig]  =  ex(8)  [s3Ix  Ig+82DM  s(K  Ig  +  Hr2)  +  K  d] 
The  characteristic  equation  of  the  system  is  the  re  Core 


-  2_  (K  I  +  H  )  K  D 

3  .  s  D  f  _g _ _ 


8  +  - 


+  8  - 


=  0 


(3.4) 

(3.5) 

(3.6) 

(3.7) 

(3.8) 

(3.9) 

(3.  10) 


Jx  *g 


The  stability  of  this  system  can  be  demonstrated  by  applying  the  Routh  stability 
criterion.  Consider  the  following  third  order  equation- 


a  s  ^  +  bs^  +  cs  +  d  =  0 

corresponding  to  this  equation,  the  Routh  array  is 
a  c 

b  d 

bc-ad 

b 

d 


(3.  11) 


If  the  coefficients  (a,  b,  c,  d)  of  equation  (3.  11)  are  positive,  a  necessary  and 
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sufficient  condition  for  stability  is  that 

b  c  -  a  d  >  0  ( H ,  I  2  > 

applying  this  condition  to  the  characteristic  equation  (3.  10)  yields: 

I  D(K  I  +H2)-I  I  KD>0  (3.  IS) 

x  '  g  r  '  x  g 

which  simplifies  to 

I  D  H  2  >  0  (3.  14) 

x  r 

which  clearly  is  always  satisfied.  Hence  stability  of  the  system  is  assured. 

Intuitively  this  is  to  be  expected. 

Since  the  characteristic  equation  (3.  10)  has  only  real  coefficients, 
the  set  of  roots  can  only  be  of  two  types,  i.e. 

1)  three  real  roots 

2)  one  real  root  and  a  pair  of  complex  conjugate  roots- 

It  is  seen  from  the  characteristic  equation  (3.  10)  that  the  sum  of  the  roots  is  given  by: 

,,T2.r3.-  °  (3.15) 

where  D  is  the  damping  coefficient 

I  is  the  moment  of  inertia  of  the  gyroscope 
g 

rl,  *  r^  represent  the  roots  of  the  equation 

Since  the  imaginary  parts  of  the  complex  conjugate  roots  cancel  on 

addition 

Re(rjH  Re  (r2)+ Re(r3)=  -  ^  (3.16) 

g 

th 

where  Re(tr  .)  represents  the  real  part  of  the  i  root. 

f  For  a  fixed  gyroscope  moment  of  inertia  (I  )  and  for  a  fixed  damping 
coefficient  (D),  the  sum  of  the  real  parts  of  the  roots  is  constant.  As  K,  the  spring 
constant,  and  Hr»  the  angular  momentum  of  the  rotor  are  varied,  the  roots  of  the 
<.  haracteristic  equation  must  migrate  such  that  the  sum  of  their  real  parts  remains 
constant.  As  some  of  the  roots  migrate  toward  the  left,  the  others  must  move  toward 
the  imaginary  axis.  This  strongly  suggests  that  optimum  damping  is  attained,  i.e. 
steady  state  is  attained  in  minimum  time,  when  the  real  part  of  all  the  roots  are  equal. 
This  condition  can  occur  in  two  ways. 

i)  a  third  order  pole  on  the  real  axis  (see  Fig.  3-2) 
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ii)  a  real  pole  and  a  pair  of  complex  conjugate  poles  all  with  the  same 
real  part,  (see  Fig*  3-3) 


Fig.  3-2  Third  Order  Pole  on  the 
Real  Axis 


Fig.  3-3  Real  Pole  and  a  Pair  of 
Complex  Poles 


Consider  case  (i)  where  a  third  order  pole  on  the  real  axis  occurs. 
For  this  case,  the  characteristic  equation  (3.  10)  must  reduce  to  the  form 

(s  +<r)3  =  0  (3.  17) 

Expanding,  and  equating  to  the  expression  for  the  characteristic  equation  (3.  10), 
yields : 


*  ?  ?  *  i  D  (K  I  ,  H  2)  K  D' 

s  +  3  <r  s  +  3  <r  s  +  <r  =  s  +y-s  +  \ — j§ — j — — '  s  +  j — 

g  *  g  x  g 


g 

Equating  coefficients  of  like  order  terms  yields: 

D 


3  <r 


r 

g 


3(r 


,  K  I  +  H 

2  =  -p-r 


*  g 


3  KD 

<T  = 


*  g 


(3.  18) 


(3.19) 


(3- 20) 
(3.21) 
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Solving  for  <r  ,  Hy  and  K  yields : 


(3  22) 


(3.23) 


(3.24) 


Clearly,  for  a  fixed  D  and  1^  these  equations  are  parametric  in  1^  .  For 

a  specified  I  and  I  ,  <r  ,  H  and  K  can  be  plotted  as  a  function  of  D  In  order 
x  g  r 

to  obtain  an  appreciation  for  the  order  of  magnitude  of  o' ,  H  and  K  ,  graphs  are 

r  2. 

drawn  as  functions  of  D  for  a  body  whose  moment  of  intertia  is  27.7  slug -ft. 

(cee  Fig.  3-4)  .  The  moment  of  inertia  of  the  gyroscope  will  be  taken  as  one  percent 
of  the  moment  of  inertia  of  the  ve hide.  (1^  =  0.  3  slug  -  ft.  ^ ) *  In  a  satellite  vehicle^it 
is  of  utmost  importance  to  keep  the  gyroscope  as  small  as  possible  so  that  more  data 
gathering  payload  can  be  included.  The  range  of  D  chosen  is  from  0.  5  x  10  3 
to  4.  0  x  10~3  ft.  -lb  j  ~~  •  These  values  are  typical  of  the  damping  attainable  in  a 

gyroscope.  It  is  seen  from  Fig.  >4  that  <r  and  Hr  are  linear  functions  of  D» 
whereas  K  varies  with  the  square  of  D  .  In  order  to  obtain  an  estimate  of  the 
settling  time  of  the  ey sten^i*  is  necessary  to  look  at  its  free  response.  The  free 
response  of  a  system,  which  has  a  third  order  pole  on  the  real  axis,  is  of  the  form 


e(t)  =  (Ko  +  Kjt  +  K2t2)  e-  fft 


(3.25) 


_  _  1  _  ^ 

It  is  clear  that  the  term  K^e  will  decay  to  e’  or  2  %of  its  maximum  value  when 


t  =  —  sec. 

O' 

-<r  t  —4 

The  term  K,  te~  will  decay  to  e*  of  its  maximum  value  when 


(3.26) 


(3.27) 


The  solution  of  this  equation  yields: 


t  =  sec. 

(T 


The  term  K^e^*  will  decay  to  e~*  of  its  maximum  value  when 

.2  -<rt  4  -6 

t  e  =  e 

ff 


(3.28) 


(3.29) 
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The  solution  of  this  equation  yields  : 

t  =  2  sec.  (3.  30) 

cr 

To  see  which  term  predominates!  it  is  necesary  to  consider  the  transfer  function  of, 
the  system  (3.9)  .  Since  the  system  has  a  third  order  pole  on  the  real  axis^the  transfer 
function  becomes 


G  (s) 

x'  ' 


(D  +  8  I  > 

.  -  ■■  -  — Jg  .  ■ 

I  I  (8  +  (T  )3 
X  g 


Td(s) 


(3.  31) 


The  free  response  of  the  system  is  given  by  its  response  to  an  impulse  torque 
(T^  (s )  =  1)  .  For  this  torque 


D  +  s  I 

Ms)  =  - 8 - r 

x  I  I  (s  +  <r)3 

X  g  '  / 


(3.  32) 


Taking  the  inverse  Laplace  transform 

0  (t)  =  i-  te’a  4  U(t)+TrHr-t2e‘<rt  U(t)  (3.  33) 

X  x  x  g 


For 


D  = 


2  x  10“ 


3  ft-lbs 


I  =0.3  slug-ft.  ^ 

K 

I  =  30  slug-ft.  ^ 


.3  -3  <3*34) 

6x(t)  =  33.  3xl0"3te‘2' 22x10  1  U(t)+ 74x  10'6t2e‘2‘ 22x10  ‘  U(t) 


^  7  5  1  A  "  ^  i 

Clearly,  the  term  33. 3x10  t e  U(t)  predominates,  hence  it  will  be  assumed 

that  the  system  damps  in  6.9  time  constants.  In  Fig.  3-5  the  settling  time  is  plotted 
as  a  function  of  the  damping  coefficient.  The  time  is  represented  as  a  percentage  of  a 
90  minute  orbit. 


It  will  be  instructive  to  sketch  the  root-loci  of  the  system  as  a  function 

of  the  damping  coefficient  (D),  the  angular  momentum  of  the  rotor  (H r ) »  the  spring 

constant  (K),  and  the  moment  of  inertia  of  the  gyroscope  (1^)-  From  Fig.  3-4  it  is 

seen  that  when  I  =  30,  I  =  0.  3,  a  third  order  pole  occurs  on  the  real  axis  at 
x  8 

0-  =  -  2.22xl0'3  when  =  109xl<T3,  K  =  49.4xl0*6  and  D 


0.002  .  In  sketching 
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Fig.  3-5  Settling  Time  as  a  Function  of  the  Damping  Coefficient  for 
Optimization  with  a  Third  Order  Pole  on  the  Real  Axis 


( 
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the  root-loci^each  of  the  parameters  will  be  in-turn  varied  while  the  remaining 
parameters  are  held  fixed  at  the  values  given  above.  Therefore,for  each  of  the  root- 
loci,  all  the  branches  will  coalesce  at  s  =  -  2.22xl0”2 


a)  Root-locus  as  a  function  of  the  damping  coefficient  (D).  (See  Fig.  3^6)  . 


The  characteristic  equation  (3.  10)  can  be  factored  into  the  form: 


where 


+  b  s)i 


P(a8Ztc) 
s(*2+h)  | 

I 


=  0 


(3.35) 


(3.36) 


b 


K  I  +  H  2 

g  t 

I  I - 

x  g 


c  = 


(3.37) 

(3.  38) 


The  roots  of  the  characteristic  equation  (3.  10)  are  therefore  seen  to  satisfy  the 
equation 


JSlLlLtfl  =r-  1  <  (3.  39) 

s(s2  +  b) 

The  poles  and  zeros  of  the  expression  on  the  left  side  of  equation  (3.  39)  are  given  by: 


P1  =  0 


P2'  P3 


z 


V 


Substituting  the  chosen  values  yields: 


p2.  p3  =  t  j  (3.86  x  10'3) 
Zj,  z2  =  t j  (1.29  x  10*3) 


(3.40) 

(3.41) 


(3-42) 

(3.43) 


The  root-locus  can  now  be  easily  constructed. 


At  D  =  0,  the  poles  of  the  system,  are  as  expected,  on  the  imaginary 
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Im  x  I03 


Fig.  3-6  Root  Locus  as  a  Function  of  the  Damping  Coefficient  for 
Optimization  with  a  Third  Order  Pole  on  the  Real  Axis 
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axis,  and  the  system  is  unstable.  Maximum  stability  occurs  when 
*-3  ft.  lbs. 


D  =  2. 0  x  10 


rad/s 


For  this  value  of  D  ,  a  third  order  pole  exists  on  the 


sec . 


real  axis;  the  values  of  H  ,  K,  and  I  were  chosen  so  that  this  would  happen.  As 

r  g 

D  is  increased  beyond  this  value*  the  free  response  of  the  system  deteriorates.  It 
should  be  noted  that  for  very  large  D,  two  of  the  roots  return  to  the  imaginary  axis. 
This  occurs  because  for  a  very  large  damping  coefficient,  relative  motion  between 
the  gyroscope  gimbal  and  the  body  cannot  occur  and  thus  the  system  behaves 
as  a  rigid  body. 

b)  Root-locus  as  a  function  of  the  angular  Momentum  of  the  rotor  (H^)  (see  Fig.  3  *7). 

For  this  case  the  characteristic  equation  can  be  factored  into  the  form: 


where 


2  (  C  8  H  2  \ 

+  a  8  +  b  «  +  d)  l  1  +-- - 2 - - -  1  =  0 

\  s  +as  +  b  s  +  d / 

(3.44) 

D 

a  =  — 

(3.45) 

g 

b  -  r 

(3.46) 

X 

1 

c  =  rr 

x  g 

(3.47) 

.  KD 

d  *  -j-y- 

*  g 

(3.48) 

expression 


For  the  chosen  values  of  D,  I^and  K,  the  poles  and  zeros  of  the 


c  s  H 

r 

— - - 2 - 

s  +  a  s  +  b  s  +  d 


(3.49) 


become 


Zj  =  0 


=  -  6.7  x  10“  J 
p2,  p3  =  t  j  (1.29  x  10*3) 


(V  '  -I) 

O-  51) 


It  is  noted  from  the  root- locus  plot  (Fig.  3-7)  that  the  damped  radian  frequency 
increases  as  increases.  This  is  to  be  expected  since  as  can  be  seen  from  the 
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Fig.  3-7  Root  Locus  as  a  Function  of  the  Angular  Momentum  of  the 
Rotor  for  Optimization  with  a  Third  Order  Pole  on  the  Real  Axis 
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characteristic  equation  (3.  10),  the  coefficient  of  the  linear  term  is  given  by 


v  H  2 
K  x  r 

r  +rr- 

X  X  g 


(3.52) 


Therefore,  affects  the  system  in  the  same  manner  as  the  spring  term.  Hence  as 
increases,  the  system  becomes  "stiffer"  and  therefore  the  damped  radian 
frequency  increases. 

c)  Root-locus  as  a  function  of  the  spring  constant  (K)  (see  Fig.  3-8). 

For  this  case  the  characteristic  equation  (3.  10)  can  be  factored  into  the 

form; 


where 


D 

r~ 

g 


(s3  +  a  s*  + 


b.)  (h  K!e  * .♦  J)  ) 

\  s3  +a  s  +b  s  / 


(3.53) 
(3  54) 


c  = 


d  = 


H 

r 

r~r 

*  g 


X 

D 

n 

*  g 


(3.55) 

(3.56) 

(3.57) 


For  the  chosen  values  of  D,  I  and  H  ,  the  poles  and  zeros  of  the 

g  r 


expression 


become 


K(c  s  4  d) 
s3  +  asZ+bs 


Z1  = 


6.67  x  10 


-3 


Pl  =  0 


3. 34  x  10'3  t  j  (1. 37  x  10‘3) 


(3.58) 

(3.59) 
(3.61) 


p2‘  p3  ! 

For  large  values  of  Kjthe  characteristic  equation  (3.  10)  reduces  to 


3  .  2D 

K  K  D 

/  .  D\ 

/  2  ,  K\ 

s  +.  r  +  .r+rT-  = 

[•  +r) 

g 

X  x  g 

\  g/ 

\  x  / 

(3.61) 


44 


Fig.  3-8  Root  Locus  as  a  Function  of  the  Spring  Constant  for  Optimization 
with  a  Third  Order  Pole  on  the  Real  Axis 


For  large  K,  the  system  will  have  a  pole  on  the  real  axis  at  s  =  -  and  will 

/v  g 

oscillate  with  a  radian  frequency  of  w  =  / ■= — •  Hence  as  K  is  increased, the  natural 

n  w  i 
V  x 

radian  frequency  increases,  and  an  undamped  oscillatory  mode  is  approached.  This 
is  to  be  expected  since  as  K  is  increased,  the  system  becomes  "stiffer"  and 
more  energy  is  stored  in  the  spring  which  will  take  longer  to  dissipate. 

d)  Root-locus  as  a  function  of  moment  of  inertia  of  the  gyroscope  (1^)  (see  Fig.  3*9] 

For  this  case  the  characteristic  equation  (3.  10)  can  be  factored  into  the 


form: 


where 


b  = 


( s 3  +b  s)  (i  +  *  i—  y. 8  +  dj\ 

\  Jg  8 (  82  +  b)  J 

(3.62) 

D 

(3.63) 

K 

5 

(3.64) 

H  2 
r 

£ 

(3.65) 

K  D 

r~ 

X 

(3.66) 

For  the  chosen  values  of  D,  Hr  and  K  the  poles  and  zeros  of  the 

expression 


1  a  s^  +  c  s  +  d 
s(s2  +  b) 


(3.67) 


become 


Pi  =  0 

p2.  P3  =  *  j  (1.29  X  10-3) 

z.,  z2  =  -  0.99  x  10*3  t  j  (0.82  x  10‘3) 


(3.68) 


(3.69) 


When  the  moment  of  inertia  of  the  gyroscope  becomes  very  large,the  characteristic 
equation  (3.  10)  reduces  to 


4fi 


Fig.  3-9  Root  Locus  as  a  Function  of  the  Moment  of  Inertia  of  the  Gyroscope 
for  Optimization  with  a  Third  Order  Pole  on  the  Real  Axis 
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The  roots  of  this  equation  are 


(3.71) 


Therefore,  as  I  increases  the  damping  of  the  system  decreases  because  the 
gyroscope  precession  rate  becomes  small  and  relative  motion  between  the  gyroscope 
gimbal  and  the  body  is  curtailed. 

Now  consider  case  (ii)  where  optimization  is  obtained  with  a  real  pole  and 
a  pair  of  complex  conjugate  poles  all  with  equal  real  parts.  For  this  case,  the 
characteristic  equation  (3.  10)  must  reduce  to  the  form: 

(s  +  <t  )  (s  +  <r  -  j  cod)  (s  +  (r  +  j  cod)  = 

(3.72) 

s  **  +  3<r  s 2  +  (<o  2+2o-2)s+<ru)  ^  =0 

'  n  '  n 


where 


(3.73) 


Clearly  tr  cannot  be  made  larger  than  co^  putting  an  upper  bound  on  the  damping 
coefficient  that  can  be  used  in  this  type  of  optimization.  Equating  expression  (3.72) 
to  the  characteristic  equation  (3.  10)  and  comparing  like  order  terms  yields: 


(3.74) 


Solving  for 


yields : 


K  I  +H2  =  w2+2  <r  2 
g  r  n 

r  i 


x  g 
2 

>  = 
n 


K  D 

r^r 


3  K 
“I 

x 


(3.75) 

(3.76) 

(3.77) 


It  is  noted  from  the  above  equations  that  three  independent  quantities  must  be  specified. 

K  and  I  will  be  fixed,  and  H^,  cr  ,  and  to  d  will  be  obtained  as  a  function  of  D.  As  in  the 

previous  case,  the  moment  of  inertia  of  the  body  plus  gyroscope  (1^)  will  be  taken  to  be 

30- slug -ft.  and  the  moment  of  inertia  of  the  gyroscope  will  be  taken  to  be  1%  of  I  or 
2  r  X 

0.  3  slug-ft.  K  will  be  chosen  for  convenience  to  be  10"  ft.  lbs.  .  Corresponding  to 

this  value  of  K,  to  2  =  10  ^  (£ii^)  and  the  upper  bound  for  D  is 
n  sec 
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to 

o 

K 

o 

LU 

CO 

X 

Q 

< 

cr 

xt 

.94 
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Fig.  3-10  Variation  of  <r  ,  and  with  the  Damping  <Z o^fficient  for 
Optimization  with  a  Real  Pole  and  a  Pair  of  Complex  Poles 


49 


-  3  ft  •  -  lb  s . 

0*9  x  10  .  •  Using  these  values,  H  ,  <r  ,  and  w  ,  are  plotted  as  a  function 

rad/sec.  r  d 

of  damping,  (see  Fig.  3  -10). 

The  free  response  of  a  system  that  has  a  pole  on  the  real  axis  and  a  pair 
of  complex  conjugate  poles,  is  given  by: 


6(t)  =  Kj  e'<rt  +  K2  e‘Tt  cos(wdt  +  $) 


{3.78) 


The  settling  time  of  this  system  is  clearly 


-.-j 


(3.79) 


In  Fig.  3^11,  the  settling  time  is  plotted  as  a  function  of  damping  coefficient  as  a 
percentage  of  a  90  minute  orbit.  Comparing  Fig.  3-11  with  Rg.  3-5  one  sees  that 
faster  settling  times  can  be  obtained  for  the  situation  employing  a  pair  of  complex 
poles.  This  is  accomplished  at  the  expense  of  permitting  the  system  to  ring  while  it 
settles . 


As  was  done  previously,  the  root-loci  of  the  system  will  be  sketched  as 

a  function  of  the  damping  coefficient  (D),  the  angular  momentum  of  the  rotor  (H^)* 

the  spring  constant  (K)  and  the  moment  of  inertia  of  the  gyroscope  (I  )*  From  Fig. 

3 "  1 0  it  is  seen  that  for  I  =  30,  I  =  0.  3,  K  =  10  the  complex  poles  and  the  real 

x  8  .4  _  3 

pole  will  have  equal  real  parts  when  D  =  5  x  10  ,  H  =3.41x10  .  It  is  also 

r  -3 

seen  from  Fig.  3-10  that  corresponding  to  these  parameter  values.  <r  =  *5.  55  x  10 

_3 

and  o)^  =  0.835  x  10  .In  sketching  the  root-loci,  ench  of  the  parameters  will  be 
in-turn  varied  while  the  remaining  parameters  are  held  fixed  at  the  values  given  above. 

a)  Root-locus  as  a  function  of  damping  coefficient  (D)  (See  Fig.  3 .*12)  . 

The  characteristic  equation  (3 » 1 0)  can  be  factored  into  the  form  of 
equation  (3.35).  Substituting  the  chosen  parameters  into  equation^ 3.  40)  and  (3.41) 
yields : 


(3.80) 

p2,  P3  =  -  j(l-28  x  10  >) 

Zj,  z2  =  t  j(0.  57  x  10“ 3) 

(3.81) 

At  D  =  0  the  poles  of  the  system  are  on  the  imaginary  axis  and  the 
system  is  unstable.  As  D  becomes  very  large,  two  of  the  roots  return  to  the 
imaginary  axis  causing  the  free  response  of  the  system  to  deteriorate.  As  in  the  case 
of  optimization  with  a  third  »*»der  pole;the  system  begins  to  behave  as  a  rigid  body. 


50 


DAMPING  COEFFICIENT  (D)  x  I04  FT-LBs/ RAD/SEC 


Fig.  3-11  Settling  Time  as  a  Function  of  the  Damping  Coefficient  for  Optimization 
with  a  Real  Pole  and  a  Pair  of  Complex  Poles 
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Fig.  3-12  Root  Locus  as  a  Function  of  the  Damping  Coefficient  l'or  Optimization 
with  a  Real  Pole  and  a  Pair  of  Complex  Poles 
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When  D  =  5x10  all  the  roots  of  the  characteristic  equation  have  oqual  real  parts. 
This  is  the  approximate  value  of  D  where  optimum  damping  occurs.  An  examination 
of  Fig.  3-10  indicates  that  the  root-locus  appears  tangent  to  cr  =**5.55  x  10"^  .  In 
actuality  however,  there  is  a  range  for  D  where  all  the  roots  of  the  characteristic 
equation  lie  to  the  left  of  the  line  <r  =  -  5.55  x  10"^  .  This  can  be  shown  in  the 
following  manner:  (s-fol)  is  substituted  for  s  in  equation  (3.  10)  thereby  increasing 
all  the  roots  of  the  characteristic  equation  by  Joj  .  The  resulting  equation  is; 


s^  +  s^  (-3fcr|  +-p)  +  s{5 o' ^ 


(3.82) 


=  0 


If  for  some  range  of  D  all  the  roots  of  this  equation  lie  in  the  left-half  plane,  then  for 
that  range  of  D,  all  the  roots  of  the  characteristic  equation  (3.  10)  must  lie  to  the  left 
of  the  line  <r  =  -5.55  x  10  ^  . 


Applying  the  Routh  test  to  equation  (3.82);it  is  seen  that  the  values  of  D 
which  place  the  roots  in  the  left-half  plane  must  satisfy  the  inequality 

2  ->  \ 


2kr\D 

T2' 

g 


+  D 


10  (T  ,  , »  3  .  -  2  ^  ^ 

— - -  +  1 2jcr|  +  dp\  to  >0 

I  n 


(5.83) 


If  the  left  side  of  the  above  expression  is  set  equal  to  zero,  and  the  resulting 
quadratic  equation  in  D  is  solved,  either 

1)  a  pair  of  complex  conjugate  roots 

2)  two  equal  positive  real  roots 

3)  two  distinct  positive  real  roots 

is  obtained.  The  first  case  indicates  that  all  of  the  roots  of  the  characteristic 

equation  never  exist  simultaneously  on  or  to  the  left  of  the  line  tr  =  -5.  55  x  10 

The  second  case  indicates  that  for  some  value  of  D  the  complex  pair  of  roots  of  the 

characteristic  equation  lie  on  the  line  cr  =  -  5-55  x  10  \  and  the  real  root  lies  either 

on  or  to  the  left  of  the  line.  The  third  case  indicates  that  there  is  some  range  of  D 

where  all  the  roots  of  the  characteristic  equation  lie  to  the  left  of  the  line 
-3 

cr  =  -  5.55  x  10  For  the  equation  under  consideration,  the  two  roots  are  indeed 

positive  and  distinct. 


D,  =  I 
1  g 


+  dm 


(3.84) 


D2  =  *1  Ig 


(3.85) 
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Fig*  3-13  Root  Locus  as  a  Function  of  the  Angular  Momentum  of  the  Rotor  for 
Optimization  with  a  Real  Pole  and  a  Pair  of  Complex  Poles 
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Hence  it  is  assured  that  there  is  a  range  of  D  where  a  better  settling  time,  than 
predicted  by  the  optimization  procedure,  is  possible.  From  Fig.  3—12  it  is  obvious  that 
the  difference  is  slight  and  the  optimization  procedure  will  give  very  close  to  the  best 
possible  settling  time.  The  reason  for  the  optimization  procedure  not  yielding  the 
smallest  settling  time  is  that  the  procedure  assumed  that  the  sum  of  the  roots  remained 
constant.  This  is  not  true  since  the  sum  of  the  roots  depends  on  D  ,  and  D  is  varying 
(see  equation  3.  16) 


b)  Root-locus  as  a  function  of  the  angular  momentum  of  the  rotor  (Hr).  See 
Fig.  3-13). 

For  this  case  the  characteristic  equation  (3.  10)  can  be  factored  into  the 
form  of  equation  (3.44).  Substituting  the  chosen  parameters,  equations  (3.^0)  and 
(3.  51)  become 


Pl  =  -  1.667  x  \0~J 

P2>  P3  =  -  j  (0- 575  X  10'3) 


(3.86) 

(3.87) 


As  H  is  increased,  th©  damped  frequency  of  oscillation  increases.  The  explanation 
r 

for  this  is  identical  to  that  given  for  the  case  of  optimization  with  a  third  order  pole 
when  the  corresponding  root  locus  was  discussed. 

c)  Root-locus  as  a  function  of  the  spring  constant  (K),  (see  Fig.  3—14). 


For  this  case  the  characteristic  equation  (3.  10)  can  be  factored  into  the 
form  of  equation  (3,53)  .  Substituting  the  chosen  parameters,  equations  (3.59)  and 
(3. 60)  becomes 

z  |  =-1.67  X  10'3  (3.88) 


Pl  =  ° 

P2’  P3  = 


(3.89) 

-  0.833  x  10*3  -  j  (0.77  x  10*3) 


The  analysis  of  the  system,  for  large  values  of  K  is  identical  to  that  described  for 
the  case  of  optimization  with  a  third  order  pole  when  the  corresponding  root-locus 
was  discussed. 

d)Root-locus  as  a  function  of  gyroscope  moment  of  inertia  (1^).  (see  Fig.  3-15) 

For  this  case,  the  characteristic  equation  (3.  10)  can  be  factored  into  the 
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form  of  equation  (3.  62)  . 
(3 . 69)  become 


Substituting  the  chosen  parameters,  equations  (3.68)  and 

Pi  =  0 

(3.90) 

p2,  p3  =  -  j  (0.577  xlO'3) 

Zj,  z2  =  -  0.  387  x  10‘3  t  j  (0.  43  x  10'3)  (3. 91) 


It  is  seen  from  Fig.  3-15  that  there  is  a  range  of  I  where  the  response  of  the  system 
is  better  than  that  obtained  by  the  optimization  procedure.  However,  the  optimization 
procedure  will  yield  a  value  which  is  reasonably  close.  The  reason  for  the  optimization 
procedure  not  yielding  the  best  result  is  that  the  sum  of  the  roots  depend  on  1^  and 
therefore,  for  a  varying  1^,  is  not  constant.  (See  equation  (3.16))  .  The  analysis  of 
the  system  for  large  values  of  1^  is  identical  to  that  described  for  the  case  of 
optimization  with  a  third  order  pole  when  the  corresponding  root-locus  was  discussed. 

C)  Analysis  of  a  Simple  Configuration  Employing  Direct  Gyroscope  Damping 

Direct  Damping  can  be  demonstrated  with  the  following  example.  Consider 
the  geometry  shown  in  Fig.  3-16.  In  the  geometry  pictured,  the  body  is  constrained 
to  rotate,  with  a  spring  type  restoring  torque,  about  the  z  axis  of  the  body.  The 
gyroscope  is  again  constrained  to  precess  relative  to  the  body  about  the  z  axis. 
Consider  a  disturbance  torque  applied  to  the  body  about  its  z  axis.  In  this  case  one 
obtains 

T  .  =  I  V  +  K  6  +  D(6  -  6  )  (3-92) 

a  z  Z  z  Z  g 

where  I  is  the  moment  of  inertia  of  the  body  about 

the  z  axis  of  the  body 
K  is  the  spring  constant 

D  is  the  coefficient  of  damping  between  the  body 
and  the  gyroscope  gimbal. 


The  torque  applied  to  the  gyroscope  is  given  by 


D(G 

z 


I 

g 


e 

g 


(3.93) 


I  is  the  moment  of  inertia  of  the  gyroscope 
g 

about  the  z  axis  of  the  body. 


where 


Fig.  3-15  Root  Locus  as  a  Function  of  the  Moment  of  Inertia  of  the  Gyroscope 
for  Optimization  with  a  Real  Pole  and  a  Pair  of  Complex  Poles. 
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Introducing  Laplace  transform  notation,  one  obtains 


and 


T,  (s)  =  I  s  e  (s)  +  K  e  (s)  +  Ds  0  (s)  -  D  s  0  (s) 


D  s  ez(s)  -  D  s  0g(s)  =  I  8  eg(s) 


Solving  these  equations  for  6^(8)  yields: 


Td(s)[D+s  ig]  =  ez(8) 


*3I  I  +s2D(I  +  I)+sIK+KD 

z  g  g  Z  g 


(3-94) 

(3.95) 

(3-96) 


(3.97) 


Therefore  the  characteristic  equation  for  a  system  employing  direct  damping  is: 

3  ,  .2  D(Iz+  V±  sK  j_  KD  n 
s  -f  s  — I  *1 —  +  I — T  *  0 

z  g  z  z  g 

Applying  the  Routh  stability  criterion  yields  as  a  necessary  and  sufficient  condition 
for  stability 


D(I  +  I  )  K  I  -  I  I  K  D  >  0 
'  g  z  g  z  g 


(3.98) 


which  reduces  to 

K  D  Ig2  >  0 

which  is  clearly  always  satisfied.  Therefore  the  system  is  always  stable.  This  was 
intuitively  expected. 

As  in  the  case  of  indirect  damping,  the  characteristic  equation  (3. 9*7) 
has  only  real  coefficients.  Therefore,  the  set  of  roots  can  only  be  of  two  types,  i  e 

1)  three  real  roots 

2)  one  real  root  and  a  pair  of  complex  conjugate  roots. 

The  sum  of  the  real  part  of  the  roots  of  equation  (3.97)  is  given  by: 

D(I  +  I  ) 

Re(r,)+ Re(r2)  Re(r3)=-  — — —  (3.100)' 

z  g 

For  a  given  body  (fixed  lz),and  for  a  fixed  D  and  1^ ,  the  sum  of  the  real  part  of 
the  roots  is  constant.  This  strongly  suggests  that  optimum  damping  can  be  obtained 
through  the  same  procedures  used  for  indirect  damping.  It  should  be  noted  that  the 
angular  momentum  of  the  rotor  (Hf)  does  not  enter  into  the  characteristic  equation 
(3.97),  hence  there  are  less  parameters  that  can  be  adjusted  to  meet  the  required 
constraints  than  for  the  case  of  indirect  damping.  It  will  be  shown  that  neither  of  the 
two  optimization  methods  will  yield  any  realizable  results. 
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Fig.  3-16  Configuration  to  Illustrate  Direct 
Gyroscope  Dampiifg 
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Optimization  with  a  third  order  pole  on  the  real  axis  will  be  considered 
first.  The  characteristic  equation  (3.97)  must  again  reduce  to  the  form  of  equation 
(3,17).  Therefore: 


a  j  2  3 

s’5  +  3(r  8  +  3a  s+a 

3  2  D(Ig  +  V  .  SK  *  KD 

=  ■  +  s  — f-j —  +-r  +  rr 

g  z  z  *  8 

(3.101) 

Equating  coefficients  yields: 

D(I  +  I  ) 

3«r  =  ;S  t  . 

g  z 

(3102) 

.2  K 

3a  =  r 

z 

(3.103) 

3  KD 

a  "  T-r 

Z  g 

(3.104) 

A  solutibn  to  the  above  equation  is 

only  possible  if 

I  =  81 

(3.105) 

g  * 

Clearly,  this  is  not  a  practical  situation;  optimization  of  this  form  is  not  attainable. 


Now  consider  optimization  with  one  real  root  and  a  pair  of  complex 
c  onjugate  roots.  The  characteristic  equation  (3.97jmust  reduce  to  the  form  of 
equat ion  (3.  72).  Therefore: 


3  .  ,  2  .  .  2  2.  .  2  3  2  D(Ig  +  V 

i  +3(T8  +(co  +  2a  )  8  +  a  W  =8  +S  - V  - - 

n  '  n  IT  I 


g 


,  sK  ,  KD 

+  t~  +r -r 

z  z  g 


Equating  like  order  coefficients  yields: 


3. _  d,i6  *  y 

•  a 

2  .  ^  2  K 

w  +  2a  -  -r- 

n  T 


2  KD 

a“n  *  TT 
z  g 


Combining  equation  (3.107)  and  (3.109)  one  obtains 

2  3K 

%  -  T-*r 

g  z 


(3.106) 


(3.107) 

(3.108) 

(3.109) 


(3. 110) 
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Substituting  equation  (3.110)  into  equation  (3.108)  yields: 


K 

r 

z 


3K 
g  Z 


=  2 <r  ' 


Since  <r  must  be  greater  than  zero. 


K 


z 


3K 

TT 


(3.111) 


(3.112) 


or 


Ig>2Iz  (3.113) 

As  for  the  case  of  optimization  with  a  third  order  pole  on  the  real  axis,  optimization 
with  a  real  pole  and  a  pair  of  complex  conjugate  poles  requires  the  moment  of  inertia 
of  the  gyroscope  to  be  larger  than  the  moment  of  inertia  of  the  vehicle  itself.  This  is. 
of  course, not  practical. 


In  order  to  obtain  some  appreciation  as  to  the  nature  of  direct  damping, 
root -loci  of  the  system  will  be  sketched  as  a  function  of  the  damping  coefficient  (D), 
the  spring  constant  (K)  and  the  moment  of  inertia  of  the  gyroscope  (1^).  The  following 
values  will  be  chosen  for  the  parameters: 


D  =  2  x  10 


ft. -lbs . 
rad/sec. 


K  »  49.  4  x  10‘6  ft.  -lbs/ rad 
1^  =  30  slug -ft? 

1^  *  0.  3  slug -ft? 


These  values  correspond  to  those  used  for  optimization  of  indirect  damping  with  a 
third  order  pole. 

a)  Root-locus  as  a  function  of  damping  coefficient  (D).  (See  Fig.  3  *17). 

For  this  case  the  characteristic  equation  (3.97)  can  be  factored  into  the 

form: 


(s 


3 


+  bs) 


P(aa2+c) 

s(s^+b) 


0 


(3.114) 
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Fig.  3-17  Root  Locue  as  A  Function  of  the  Damping  Coefficient  for 

Direct  Damping 


where 


a 


(3.11S) 
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I  +  I 
-  g  z 

“  I — I 

g  z 


b  = 


K 

I” 

z 

K 


z  g 

Substituting  the  chosen  values ,  the  poles  and  zeros  of  the  expression 

2  . 

as  +  c 


+b) 


become 


Bj,  z2  =  -  j  (1.28  x  10'3) 


Pi  *  0 


P2.  P3  =  t  j  (1.29  x  10' J) 


(3.116) 

(3.117) 


(3.118) 


(3  119) 

(3.120) 


It  is  seen  from  1  ig.  3“  17  that  very  little  damping  is  attainable.  When  indirect 
damping  was  used,  it  was  possible  to  obtain  considerably  more  damping  for  the  same 
parameters.  The  reason  for  this  is  that  in  the  indirect  case,  relative  motion  between 
the  gyroscope  gimbal  and  the  body  was  produced  as  a  result  of  gyroscopic  precession. 
The  velocity  gain  of  the  gyroscope,  the  ratio  of  precession  velocity  to  input  velocity, 
is  available  to  give  a  large  relative  velocity  for  a  small  input  velocity*  With  direct 
damping  however,  relative  velocity  between  the  gyroscope  ginbal  and  the  body  is 
the  same  as  the  input  velocity,  hence  no  velocity  gain  is  present. 

b)  Root-locus  as  a  function  of  the  spring  constant  (K)  (See  Fig.  3*18) 

For  this  case  the  characteristic  equation  (3.97)  can  be  factored  into 

the  form: 


(b3  +  a.2)  (\  +K(b8  tSl\  (3.121) 

\  s2(s+a )/ 

where 


+  u 

g  * 

(3.122) 

I  I 

g  z 

1 

T* 

(3.123) 

D 

(3.124) 

I  I 

z  g 

c  = 


y  >(D 

-6.74  -6.67 


Re  x 


Fig.  3-18  Root  Locus  as  a  Function  of  the  Spring  Constant  for  Direct  Damping 
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Substituting  the  chosen  values,  the  poles  and  zeros  of  the  expression 

(bs  +  c) 
s^(s  +  a) 

become 

=  -  6.67  X 10"3 

Pi'  P2  =  0 

p3  =  -  6.74  x  10'3 


(3.12^) 


(3.126) 


(3.127) 


c)  Root-locus  as  a  function  of  the  moment  of  inertia  of  the  gyroscope  (Ig) 

(See  Fig.  3-19). 

For  this  case  the  characteristic  equation  (3.97)  can  be  factored  into  the 

form: 


(8  +  as  +  cs)  l  1  +  -8- 


■j-  (bs2  +  d) 


b(b*  +  as  +  cj 

Substituting  the  chosen  values,  the  poles  and  zeros  of  the  expression 

(bs2  +  d) 

2 


(3.128) 


(3.129) 


8  ( s  +  as  +  c) 


become 

zv  z2  =  t  j  (1.29  X  10'3) 

p1=  0 

p2,  p3  =  t  j  (1.  28  x  10"3) 


(3.130) 

(3.131) 


Comparing  indirect  damping  with  direct  damping,  it  is  seen  that  a  much 
higher  damping  coefficient  is  required  to  obtain  a  reasonable  settling  time  for  the 
directly  damped  case-  In  a  satellite  vehicle  employing  a  single  gyroscope,  whose 
precession  axis  is  along  one  of  the  body  axes  of  the  satellite,  direct  damping  will 
occur  on  one  of  the  axes  whereas  indirect  damping  will  occur  on  another.  If  the 
damping  coefficient  is  chosen  to  provide  adequate  damping  for  the  directly  damped 
axis,  the  damping  on  the  indirectly  damped  axis  will  be  well  beyond  the  optimum 
point,  hence  poor  settling  time  about  that  axis  will  be  obtained.  This  suggests  that 
either  the  gyroscope  precession  axis  be  offset  from  the  body  axes  or  several  gyro¬ 
scopes  be  employed  so  that  indirect  damping  can  be  provided  about  each  of  the  axes 
to  be  controlled. 
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In  order  to  analyze  the  dynamics  of  a  satellite  containing  internally 
mounted  gyroscopes,  it  will  be  first  necessary  to  analyze  the  dynamics  of  a 
gyroscope  mounted,  in  an  arbitrary  position,  within  an  orbiting  vehicle.  This  will 
be  accomplished  in  Chapter  IV. 


Fig.  3-19  Root  Locus  as  a  Function  of  the  Moment  of  Inertia  of  the  Gyroscope  for 

Direct  Damping 
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CHAPTER  IV  -  GYROSCOPE  DYNAMICS 


It  was  noted  in  the  examples  at  the  end  of  Chapter  II  that  an  orbiting 
vehicle  would,  at  best,  have  only  an  oscillatory  free  response.  Any  practical  v  hide 
must,  of  course,  have  damping  which  reduces  disturbance  oscillations  in  a  reasonable 
period  of  time.  Since  a  conventional  form  of  damping,  atmospheric  drag,  present 
in  vehicles  which  travel  in  the  atmosphere,  is  not  available,  some  other  means  of 
damping  must  be  provided.  It  was  noted  in  Chapter  III  that  a  gyroscope  with  a 
damping  fluid  constraint  may  be  used  to  provide  this  damping.  In  this  chapter  the 
equations  of  motion  of  a  gyroscope  mounted  within  an  orbiting  vehicle  will  be  developed 
for  an  arbitrary  orientation  of  the  precession  and  spin  axes  relative  to  the  coordinate 
system  of  the  vehicle.  From  this  general  derrvatbn  special  cases  will  be  cited. 


A)  The  Gyroscope  Equations  for  Symmetrical  Design 

The  equations  of  motion  of  a  gyroscope  are  developed  in  Appendix  B  of 
this  report.  They  are  repeated  here  for  convenience. 

T  =  I  i)  +  (I  -  I  )  co  a)  -  H  «  (B  9) 
gi  Xc  gi  Zs  ye  g‘  gk  r  gk 
«  g  g  g  e  *g  g  g 


(B- 10) 


T  =  H  +  I  ci  +  (I  -  I  )  u  w- 

gJ  r  yg  gL  Xg  Zg  gi  gk 

g  g  g  g 

T  =Iu.  +(I>I)uu+Hw  (B.  11) 

gk  2g  8k  yg  xg  % 

g  g  g  g  g 


where  X^,  Y^»  are  the  principal  axes  of  the  gyroscope  (sp(*  Fig  Bd) 

T  »  T  ,  T  are  the  components  of  torque  on  the  gyroscope  along 
g  i  S|  g 

g  Jg  g  the  X,  Y,  and  Z  axes  of  the  gyroscope  respec  tively 


I  ,  I  ,  1^  are  the  polar  moments  of  inertia  of  the  gyros  ope  along  the 

g  g  g  X,  Y  and  Z  axes  of  the  gyroscope 

w  ,  w  cj  are  the  components  of  the  angular  velocity  of  the  gimbal 

Ri  gi  \  along  the  A,  Y  and  Z  axes  of  the  gyroscope  with  respect  to 

g  Jg  g  inertial  space 

is  the  angular  momentum  of  the  rotor 

These  equations  are  sim^ar  to  Euler’s  equations  of  motion  except  that  the 
angular  momentum  of  the  gyroscope  rotor  introduces  some  additional  terms  The 
angular  momentum  of  the  rotor  will  normally  be  held  constant  through  the  i^e  of  an 
appropriate  speed  control  system.  Therefore  can  be  taken  to  be  zero.  Th^ 

equations  reduce  to  a  particularly  simple  form  if  symmetrical  design  of  the  gyroscope 

is  assumed.  (I  =1  =1  =  I  J.  In  order  to  facilitate  the  analysis,  symmetrically 

xg  yg  zg  8 
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designed  gyroscopes  will  be  used  throughout  this  report.  It  is  anticipated  that  this 
restriction  will  not  cause  a  significant  loss  in  generality.  Making  the  appropriate 
substitutions,  equations  (B  9),  (B.  10)  and  (B-ll)  become 


T 

I 

co  -  H 

CO 

(4.1) 

g 

g.  r 

6i 

8k 

g 

g 

g 

T  = 

I 

g 

CO 

g. 

(4.2) 

Jg 

Jg 

T  = 

I 

u  +  H 

CO 

(4.3) 

gk 

g 

gk  r 

*i 

g 

g 

g 

It  is  further  assumed  that  the  body,  to  which  the  gyroscope  is  affixed,  is  reasonably 
aligned  with  the  reference  coordinate  system  so  that  the  equations  developed  for 
small  angular  deviations  in  Chapter  II  can  be  applied. 


B)  The  Transformation  Matrix  Relating  Body  and  Gyroscope  Coodinates 

In  the  derivation  to  follow  in  Section  C»the  gyroscope  will  be  assumed 

mounted  at  an  arbitrary  position  within  the  vehicle.  The  coordinate  system  of  the 

gyroscope  will  be  denoted  by  i^,  and  k^  where  i^  is  along  the  precession  axis 

of  the  gyroscope,  j  is  along  the  angular  momentum  of  the  rotor  and  k  is  chosen 
g  8 

to  forma  right  handed  coordinate  system.  (See  Fig.  B-l).  The  gyroscope  coordinate 

axes  can  be  related  to  the  body  coordinate  axes  through  an  appropriate  array  of 
Eulerian  angles.  This  procedure  is  similar  to  that  used  in  Chapter  II,  (Section  C)» 
where  a  transformation  matrix,  relating  the  body  coordinate  axes  to  the  reference 
coordinate  axes  was  obtained.  Though  the  transformation  matrix,  (2.  3)  can  be 
applied  to  the  present  situation,  it  will  not  represent  an  intelligent  choice.  It  is 
advantageous  to  choose  the  set  of  Eulerian  angles  so  that  the  location  of  the  precession 
axis  fixes  two  of  the  Eulerian  angles,  and  precession  of  the  gyroscope  causes  only  the 
third  Eulerian  angle  to  vary.  Examination  of  equation  (2.  3)  reveals  that  if  the  trans¬ 
formation  matrix,  developed  to  relate  the  body  coordinate  axes  to  the  reference  coor¬ 
dinate  axes,  is  used  to  relate  the  gyroscope  coordinate  axes  to  the  body  coordinate 
axes,  precession  of  the  gyroscope  will  in  general  cause  all  three  Eulerian  angles  to 
change.  Precession  of  the  gyroscope  represents  a  rotation  about  the  i^  axis  of  the 
gyroscope;  it  would  therefore  be  wise  to  perform  this  rotation  last.  The  Eulerian 
angles  that  will  be  used  will  be  denoted  by  <|>  ,  (f)^  and  4>^  .  These  angles  are  defined 
by  successive  rotations  of  the  gyroscope  axes  in  the  order  4>z»  4>x  which  rotate 

the  gyroscope  axes  from  coincidence  with  the  body  axes  into  its  actual  position.  (See 
Fig*  4-1).  These  angles  will  be  taken  as  positive  when  the  sense  of  rotation  is 
equivalent  to  that  of  a  right  handed  screw  advancing  in  the  direction  of  the  axis  of 
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rotation.  Fig.  4-1  shows  the  three  successive  rotations  and  their  corresponding 
transformation  matrices.  The  overall  transformation  matrix  relating  the  components 
of  a  vector  in  body  coordinates  to  the  components  of  the  same  vector  in  gyroscope 
coordinates  is  given  by 


(4.4) 


where 


A  = 


(A)  (A)  (A  ) 

x  z  y 


(4.5) 


and 


A  = 


COS  COS 


sm  <f>^  cos  <J)z 


cos  d>  sin  <b  sm  d> 
x  ^y  z 

+  sin  <b  cos  <b 
x  y 


-  sm  d>  sin  <b  sin  <b 
x  y  z 

+  cos  <b  cos  <b 

*  y 


(4.6) 


-  cos  <b  cos  <b  sin  A  ' 

x  ^y  1  z 

+  sin  6  8 in  <b 

^x  xy 


sin  <b  cos  A  sin  <b 
Xx  xy  z 

+  cos  cb  sin  <b 
x  y 


sin  4) 


sin  A  cos  A 
y  z 


cos  6  cos  A 
y  Yz 


As  for  the  transformation  matrix  relating  the  body  and  reference  coordinates,  matrix 
multiplication  is  not  commutative  and  must  be  performed  in  the  order  indicated.  This 
again  signifies  that  the  order  in  which  the  rotations  are  performed  is  of  utmost 
importance.  It  should  be  noted  that  the  first  two  rotations,  through  the  angles 
and  <j>z  ,  determine  the  location  of  the  precession  axis  of  the  gyroscope.  Precession 
of  the  gyroscope,  merely  changes  the  remaining  Eulerian  angle  <t>x  . 

C)  The  General  Gyroscope  Equations 

The  angular  velocity  of  the  gyroscope  with  respect  to  inertial  space  is  the 
vector  sum  of  the  angular  velocity  of  the  gyroscope  with  respect  to  the  body  and  the 
angular  velocity  of  the  body  with  respect  to  inertial  space. 

u>  /Inertial  =  to  /  ,  ,  +  to,  ,  /Inertial  (4.7) 

g/  space  g /  body  body/  space  '  7 

The  angiiar  velocity  of  the  body  with  respect  to  inertial  space,  in  terms  of  body 
coordinates,  is  given  as 

“b/i  =  wb.  +v  4  +  x  kb 

'  lb  Jb  kb 


(4.8) 
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where 


“b.  = 

c  +  f;  ft 

x  z  y 

(4.9) 

COi  =* 

G  +  ft 

(4.10) 

b . 

y  y 

V 

l  -  on 

z  x  y 

(4-11) 

These  equations  were  obtained  from  equations  (2.14),  (2.15)  and  (2.16)  by  setting 

Q  =  25  0  -  Using  the  vector  transformation  matrix,  derived  in  the  previous  section, 

x  z 

(equation  4.6),  one  can  express  the  components  of  the  angular  velocity  of  the  body 
with  respect  to  inertial  space  in  terms  of  gyroscope  coordinates.  The  resulting 
equation  then  becomes 


where 


»„/  -  «...  ‘g+“b 


>/.  * 

'  1  R 


bi 


4,  (® ,  +  «v>  -  Sin 


x  y  y 


cos  4>z  (6Z  -  «y> 


+  COS  <t»y  COS  <t>z  (6X  +  e7.  V 


(4.12) 


(4.13) 


s  4  cos  4  (e  +  ft  )  +  (cos  Sin  4  sin  4Z 


z  y 


J8  +  sin  *  cos  4,  )  (6  -  C  ft  )  +  (-  cos  4X  cos  4  sin  4 


x  y 

+  sin  4>x  sin  <J>y)  (f'x  +  czny) 


(4.14) 


=  -  sin  4>x  cos  4>z  (6  +  fty)  +  (cos  4>x  cos  4y 

-  sin  <(»x  sin  <)>y  sin  (6g  -  6x  fty)  +  (cos  ^  sin  4. 

+  sin  4  cos  6  sin  4)  (C'x  +  ez  B  ) 
x  y  ^  / 


(4.13) 


Since  the  gyroscope  is  constrained  to  rotate,  relative  to  the  vehicle,  about  its  pre¬ 
cession  axis  (ig)  ,  the  ang 
body  can  be  expressed  as 


cession  axis  (ig)  ,  the  angular  velocity  of  the  gyroscope  gimbal  with  respect  to  the 


=  $  i 


x  g 


(4.16) 
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Therefore,  from  equation  (4.7),  the  angular  velocity  of  the  gyroscope  gimbal  with 
respect  to  inertial  space  is  given  by 


w  /i 

g/ 


space 


where 


gi 

8 

i  +  w  j  +o>  k 

g  g:  g  gk  g 

Jg  g 

(4  17) 

u> 

gi 

g 

=  “b.  +  *x 

*g 

(4.18) 

(0 

gi 

=  “b.  . 

(4.19) 

Jg 

Jg 

0) 

gk 

=  UK 

bk 

(4.  20) 
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Differentiating  the  components  of  equation  (4.17)  with  respect  to  time,  yields  the 

components  of  the  angular  acceleration  of  the  gyroscope  gimbal  with  respect  to 

inertial  space  in  terms  of  gyroscope  coordinates.  It  should  be  noted  that  the  location 

of  the  precession  axis  is  fixed  for  a  given  gyroscope  orientation,  therefore  Ay  and 

A  are  fixed  constants  and  the  only  Euler ian  angle  that  will  vary  in  time  is  A 
z  x 

(This  angle  will  vary  as  the  gyroscope  precesses). 


*  sin  4>  (©  +  0  )  -  sin  A  cos  A(0-GD-Gf2) 

z  y  y  Yy  Yz  z  x  y  x  y' 

+  cos  4>  c  08  A  (0  +  0  12  +  0  12  )  +  A 

y  z  x  z  y  z  y  x 
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=  -  <t>  sin  *  cos  <t>  (e  +  n  )  +  cos  +  cos  A  (6  +  fi  ) 

x  x  z  y  y  x  z  y  yf 
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'  Yx  Yx  Yy  Yx  Yx  YyM  z  x  y' 

+  (  cos  A  sin  A  sin  A  +  sin  A  cos  A  )  (0  -  0  12  -  6  12  ) 
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+  (A  sin  A  cos  A  sin  A  +  A  cos  A  .  sin  A  )(G  +  G  12  ) 
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'  YX  Yx  y  Yx  X  y  Z  z  x  y 
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Substituting  equations  (4.18),  (4.19),  (4.20),  (4.21),  (4.  22)  and  (4.  23)  into  equations 
(4.1),  (4.2),  and  (4.3)  yields  for  the  components  of  torque  on  the  gyroscope  gimbal 
along  the  gyroscope  axes 


ig  [-  *x  cos  cos  4>z  (ey  +  ny)  -  Sin  <f,x  cos  +,(0+6) 

-  (  <£x  sin  <()x  COS  <b  +  <£  cos  *  Sin  A  sin  $  )(C  -  0  fi  ) 
y  x  y  z  z  x  y 

+  (  cos  4>x  cos  <f,y  -  sin  <J>x  sin  <j)y  sin  &x«y  -  fi  ) 

+  <  *  8ln  'f’x  8in  +  *x  cos  ‘•’x  cos  <*\,  sin  ♦  )(0  +0  U  ) 

/  A  y  Z  X  Z  y 

+  (  cos  d>x  Sin  6  +  sin  A  cos  A  sin  4»  )( 0  +  6  fi  +  0  ft  )1 

IP  ..  r  *  r  z  *  z  y  7.  y'j 

+  Hr  [sin  4>  (0  +0  )  -  sin  A  cos  A  (0  -  6  O  ) 

y  y  y  z  z  xy' 

+  cos  A  cos  A  (6  +  0  n  )+  A  ] 
y  z,  x  z  y'  TxJ 


(4.28) 


It  is  now  necessary  to  express  the  component:  of  torque  on  the  gyroscope 
gimbal  in  terms  of  the  body  coordinate  axes.  This  is  accomplished  by  utilizing  the 
inverse  transformation  of  equation  (4.6) 
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cos  <t>  cos  4> 

X  Z 
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Utilizing  the  inverse  matrix. 
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(4. 27) 


(4.28) 


expanding  and  performing  algebraic  simplifications,  one  obtains  for  the  components 
of  torque  on  the  gyroscope  gimbal  along  the  body  coordinate  axes 


T  =l[-4>sin4>  cos  4>  ( G  +  17  )  +  ( 0  +  0  Q  +  6  17  ) 
g  (  g*  'x  Ty  'z'y  y'  '  x  zy  z  y  ' 

^  +  (  cos  4>  cos  4>  4>  )  -  4>  sin  4>  (0  -0R)| 

y  z  Yx'  Yx  Yz  z  x  y' J 

+  H  (sin  4>  cos  4>  +  cos  4>  sin  4>  sin  6  )(0  +  f7  ) 

r  L  x  y  x  y  YzM  y  y' 

-  cos  4>  cos  4>  (6  -  G  17  )  +  (cos  4>  sin  4>  + 

x  Yz  z  x  y'  '  Yx  Yy 

sin  4>x  cos  4>y  sin  4>z)  4>x] 


(4.  2'*' 


I  T(0  +  U  )  +  sin  4>  4>  +  4>  cos  4>  cos  4>  (0 

g  l'  y  y'  Yz  Yx  Yx  y  z  z 


g  n  ) 

x  y 


+  A  sin  4>  cos  4>  (6  +0  17  )  +  H  I  (sin  4>  sin  4> 

Yx  Yv  Yz  '  x  z  y'  J  r  ^  y 


X 

-  cos 


t>x  cos  4>y  sin  4>Z)(0Z  -  0X  Hy)  -  (sin  4>x  cos  4>y 
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T  =  I  f  -  i  cos  (h  cos  6  (C  +  S1  )  +  (0  -On  -  f  U  ) 
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The  equacion  of  constraint  relating  the  gyroscope  precession  angle  (6^) 
to  the  body  angular  deviations  (0^,  and  0^)  can  be  obtained  from  the  fact  that  the 

component  of  torque  on  the  gyroscope  gimbal  along  the  X  axis  of  the  gyroscope  is 
given  by 

T  =  -  DMT  (4.  32) 

gj  x  c 

g 

where  D  represents  the  damping  coefficient 

Tc  represents  a  control  torque  applied  to  the  gyroscope. 


Equating  the  right  side  of  equation  (4.  32)  to  the  component  of  torque  on  the  gyroscope 
along  the  i^  direction  (equation  4.24)  yields  for  the  equation  of  constraint 

(4.33) 
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Equations  (4.29),  (4.  30),  (4.  31)  and  (4.  33)  are  completely  general.  For 
appropriate  values  of  6^  and  the  equations  for  any  desired  orientation  of  the 
precession  axis  can  be  obtained.  In  Section  D,  the  equations  for  several  different 
gyroscope  configurations  will  be  tabulatcd*though  by  no  means  will  this  list  be 
exhaustive . 


D)  The  Gyroscope  Equations  for  Particular  Configurations 

Although  the  gyroscope  equations  (4.29),  (4.  30),  (4.  31)  and  (4.  33)  are 
completely  general, it  will  be  convenient  to  reduce  the  gyroscope  equations  for  some¬ 
what  more  restrictive  configurations.  Several  of  the  more  important  configurations 
are  listed  below. 
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1)  The  gyroscope  precession  axis  is  located  in  the  Xb  plane,  (see  Fig.  4>-2) 

For  this  configuration,  <j>^  =  0  and  the  location  of  the  precession  axis 
is  determined  by  .  <|>  ,  of  course,  determines  the  location  of  the  spin  axis. 

Setting  $  =  0°  in  equations  (4-29),  (4.  30)  and  (4.  31)  yields  the  equations  for  the 

components  of  torque  on  the  gyroscope  when  the  precession  axis  is  located  arbitrarily 
in  the  Xb  plane 


T  =  I  [  0  +  +0fi  +  cos  <)>  if  -  (j>  sin  <b  (0  -  0  )  1  (4.  34) 

g.  gL  x  z  y  z  y  Tz  x  Yx  Yz  z  x  y'  J 

ip*  * 

b  +H  sin  <b  (G  4~  Q  )  -  cos  d>  cos  <b  (0  -  G  Q  ) 

r  [  x  y  y  x  ^z  z  x  y 

+  sin  cb  sin  A  <b  I 
x  Tz  ^x  J 


T  =  I  [ 0  +  S7  +  sin  <fi  <b  4-  6  cos  <b  (G  -  0  Q  )1 

g.  gly  y  z  x  x  z  z  xy'J 


+  H 


[  -  cos  cb  sin  (j>  (G  -  G  Q  )  -  sin  cb  (0  +0  U  ) 

r  t  rx  z  z  x  y  x  x  z  y 


-  <b  sin  <b  cos  4> 
x  z 


(4.35) 


«k. 


i  f  -  <j>  cos  a  ( e  +n)  +  e  -  en  -  e  n 
gL  *  y  y'  z  x  y  x  y 

+  i  sin  d)  (6  +  e  n  )  1  +  H  [ cos  <b  cos  <b  (6  +  t  n  ) 

Yx  ^  z  v  x  z  y'  J  r  L  Yx  z  '  x  z  y 


+  cos  <b  d>  4-  cos  cb  sin  cb  (0  +  Q  ) 

x  x  x  z  v  y  J 


(4.36) 


Setting  =  0  in  equation  (4.  3  3)  yields  the  equation  of  constraint. 


-  D  <f)  +  T  =  I  |  sin  <b  (6  4  Q  ) 

X  c  g  L  y  y' 

4-  cos  4  (e  +0  ft  4-  G  Q  )  4  A  ] 

Z  X  z  y  2  y  T  x  J 

4  H  [  sin  4>  coscb  (0  4  ft  )  -  cos  A  (0  -0  ft  ) 

r  L  ^x  Yz  y  y'  Yx  '  z  x  y 

-  sin  A  sin  (j)  (0  4-0  ft  )  1 

^X  Z  X  Z  y  J 


(4-37) 
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Fig*  4-2  Gyroscope  precession  axis  in  the  plane 


Fig*  4-3  Gyroscope  precession  axis  in  the  plane 
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2)  The  gyroscope  precession  axis  is  located  in  the  Z^  plane,  (see  Fig.  4^3) 

o 

For  this  configuration  <{>  =  0  and  the  location  of  the  precession  axis  is 

determined  by  $  •  Setting  <j>  =  0  in  equations  (4.29),  (4.  30)  and  (4.  31)  yields 

y  z 

the  equations  for  the  components  of  torque  on  the  gyroscope  when  the  precession  axis 
is  located  arbitrarily  in  the  Z^  plane. 

T  =i  \-i  sin  4>  (e  +n)  +  e  +f  ft  +6  ft  4  cos  <t>  "4  1  (4. 3«) 
g.  g  L  X  Ty  y  y'  x  z  y  z  y  7  xj 

^  4  H  i  in  <J>  cos  4>  (0  4ft)-  cos  d>  (0  -  0  ft  )  4  cos  <b  sin<(>  <b  ] 

r  L  x  y  y  y  x  z  x  y'  ^x  Ty  xj 


T  =1  fe  4ft  4  A  cos  &  (0  -  G  ft  )  +  <i»  sin  4>  ( €  4  G  ft  )1 

g  g  L  y  y  x  Yy  '  z  x  y'  x  Yy  x  z  y'J 

k  +  H  [  sin  <j>  sin  A  t  G  -  0  ft  )  -  sin  A  cos  <J>  (6  4  0ft) 

r  L  oc  ^y  z  x  y'  ^x  y  x  z  y' 


(4.  39) 


-  Sin  '•’x  ] 


T  =  I 
gk  8 
Kb 


f  -  d>  cos  A  (0  4ft  )  4  0  -  6  ft  -  8  ft  -  sin  <|)  i  ] 

l  x  y  y  y  z  x  y  x  y  y  x  j 


(4.40) 


.  .'■'V 

4  H  f  cos  <))  (0  4  0  ft  )  4  cos  <|>  cos  4>  4>  -  sin  <t>  sin  A  (6  4ft  )] 

r[  x  x  z  y  x  y  x  x  y  y  y  J 


Setting  4>z  =  0°  in  equation  (4.  33)  yields  the  equation  of  constraint 

-  D  <f>  +  T  =  I  f-  sin  4>  (0  -  0  ft  -  0  ft  )  .  4n 

x  c  g  L  y  z  x  y  x  y'  <4  41) 

+  cos  (j)  (0  4  0  ft  f  6  ft  )  4  £  1  4  H  [  sin  4>  (0  4ft) 

Yy'x  zy  zy  xj  r  L  x  y  y 


-  cos  4>  cosA(0  -  8  ft  )  -  cos  <J>  sin  <J>  (0  +6  ft  )1 

x  z  x  y'  x  y  x  z  y'J 

3)  The  gyroscope  precession  axis  is  located  in  the  Y^  plane.  (See  Fig.  4  -4) 


For  this  configuration  4>y  =  90  and  the  location  of  the  precession  axis 
is  determined  by  <(>z  .  Setting  4>y  =  90°  in  equations  (4.29),  (4.  30)  and  (4.  31)  yields 
the  equations  for  the  components  of  torque  on  the  gyroscope  when  the  precession  axis 
is  located  arbitrarily  in  the  Z^  plane. 

r  .  •  ......  .  (4.42) 

T  a  I  -  4>  cos  A  (0  4ft  )4  0  48  ft  4  0  ft  -  <|>  sin  d>  (0-0  ft)f 

g.  g  L  x  ^z  y  y'  x  z  y  z  y  Yx  z  x  y'  J 

Jb  f  *  * 

4  H  cos  <b  sin  A  (0  4ft)-  cos  A  cos  ib  (0  -  0  ft)  +  cos  A  4>  1 

r  L  oc  Tz  '  y  y'  Yx  ^z  z  x  y  x  xlj 

T  =  I  f  0  4ft  sin  <b  A  4  <?  cos  A  (0  4  0  ft  )1 4  H  [  sin  A  ( 6  -  C  ft  ) 

g,  gl  y  y  z  x  Yx  Yz  '  x  z  y'J  r  L  ^x  z  x  y' 


•  *  y 

cos  A  sin  A  (0  40  ft  )  -  d>  sin  A  cos  A 

x  z  x  z  y  x  rx  z  J 


(4.43) 
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Fig*  4-4  Gyroscope  precession  axis  in  the  plane 


BODY 


Fig*  4-5  Gyroscope  precession  axis  along  the  axis 
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T  =  I 
gk  8 

b  + 


Setting  <J)^  =  90°  in  equation  (4,3  3)  yields  the  equation  of  constraint. 


t*. 

M 


e  n  -  e  n  -  cos  <(>  $  +  *  sin*  (e  +  e  n  )1  <4,44^ 

x  y  x  y  z  x  Yz  x  z  y  J* 

cos  cf>  cos  4>  (6  +  6  R  )  -  sin  <)>  sin  cf)  -sin<j>  (0  +  R  )1) 

x  z  x  z  y'  Yx  Kz  x  Yx'  y  y  J 


,  H.45) 


-  D  4>  +  T  =  I  [sin  <f>  (0  +  R  )  -  cos  A  (0  -  G  R  .  e  R  ) 

X  c  g  L  z  y  y  Z  '  z  x  y  X  y 

+  +’  1  +  H  [sin  A  cos  A  (0  +  R  )  +  sin  A  sin  <)>  (0  -  0  R  ) 

x  J  rl  Yx  Tz  y  y'  Yx  ^z  '  z  x  y7 

-  cos  A  (0  +  G  R  )] 

Yx  x  z  y  J 


4)  The  gyroscope  precession  axis  is  located  along  X^.  (See  Fig.  4*5) 

Although  the  set  of  equations  for  this  configuration  can  be  obtained  from 

the  general  equations  (4.29)»  (4.  30)  and  (4.  31)  by  setting  <jy  =  4>z  =  0°,  they  can  also 

be  obtained  as  a  special  case  of  (1)  by  sdtting  <}>  =  tf^since,  in  that  case.  <f>  was 

z  y 

already  taken  to  be  zero.  At  any  rate,  the  resulting  equations  are 

t  =  i  ]e  +  e  n  +  e  n  +  *'  1  +  h  [  -  cos  <j>  (e  -  e  :i  )  <4-4^) 

g.  glx  z  y  z  y  -Xj  r  [  vx  z  x  y' 

°  +  sin  t  (e  +n  )] 

Tx  y  y'  J 


(4.  47) 


t  *  i  f?  +n  ♦♦  (e  -  e  c  )i  -h  rsin  <t>  (e  +  R  Vl 

gj  g[  y  y  x  z  xy  r  *  *  z  y  *  J 

Jb 

r.  (4.48) 

V  ■  *,  [-  *x  <\  ♦  V  *  ®.  -  e«  \  -  *x°y]*  r  [“«  +x  <•» *  ezny  *  ♦,»( 


% 


Setting  <j>^  =  ^2,  =  0  *n  equation  (4.  33)  or  setting  4>z  =  0°  in  equation  (4.37)  ,  yields 
the  equation  of  constraint. 


-  D  4>  +  T  =  I  f  r  +  G  R  +  0  R  +  "A  I 

x  c  v  I  x  zy  z  y  Yx  J 

+  H  [  sin  4>  (G  +  12  )  -  cos  A  (0  -  0  R  )1 

r  l  ^x  '  y  y  x  '  z  x  y'J 


(4.49) 


5)  The  precession  axis  of  the  gyroscope  is  along  the  Yb  axis.  (See  Fig.  4-6) 
As  in  (4),  the  set  of  equations  for  this  configuration  can  be  obtained  by 
setting  4>y  =  0  and  <J>^  =  90°  in  the  general  equations  (4.29),  (4.30),  and  (4.  31)  or 
as  a  special  case  of  (1)  by  setting  a  90°.  Either  way,  the  resulting  equations 
are 


gi. 


--ije 

8l  * 


e  r 

z  y 


0  R 
z  y 


(ez- 


e  n  ) 

x  y 


+  H 


[sin4>x(ey  + 


(4.-50) 


n  + 


.>] 
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:•  “^(VV'M  -Hr[co8<J,x<®z-  9xV 

Ju  •  1 


■f  sin  <)>  (0  +0  ft  )  1 

Yx  x  z  v  J 


Tgk  "  V  z 
b  + 


r  e  .  e  n  -  e  n  +  V  (e  +  e  n  )1 

L  z  x  y  X  y  Yx'x  zyJ 

H  f cos  (j)  ($  +  e  +n)l 
r  L  x  x  y  y  J 


(4.  51) 


(4. 52) 


Setting  $  =  0°  and  4>z  =  90°  in  equation  (4.  33)  or  setting  <t>z  =  90  °  in  equation 


(4.  37)  yields  the  equation  of  constraint. 


D  d>  +  T  =  I  f  e  +  h  til 
vx  c  g  [  y  y  x  J 

Hr  [CO®  K  (®z  ’  °x  V  +  ®in  ^x  (®x  +  °z  V  J 


(4.53) 


6)  The  gyroscope  precession  axis  is  along  the  axis.  (See  Fig.  4-7) 

The  set  of  equations  for  this  configuration  can  be  obtained  from  the 
general  equations  (4.29),  (4.30)  and  (4.31)  by  setting  4>y  =  -  90°  and  <t>z  =  <f, 
or  as  a  special  case  of  (2.)  by  setting  4*y  =  ”90  since  <|>z  is  already  zero.  At  any 
rate,  the  resulting  equations  are 


=  I 


-Hr[co.  *X(V  6*1!yt*x>l 


(4.54) 


i  [e  +  q  4  (nen  ) 

g  L  y  y  x  Z  y' 


lie  -  e  a  -  e  n  + 
g  II  z  x  y  *  y 

+  sin  4>  (6  +  U  )  '1 

Yx  y  y  J 


)-  Hr  tin 't’x  (®z  *  ex°y  +*x)l  (4>55) 
♦xj  +  Hr  [cos  <t>x  (ex  +  ®7  »y)  (4-  56> 


Setting  4>  =  -  90°  and  4>  =  0°  in  equation  (4.  33)  or  setting  4> 

y  z  • 

(4.41)  yields  the  equation  of  constraint. 


-D*x+Tc-  °x6y-  ®*°y+'+*S 

+  Hr  [sin  <|>x  (ey  +  ny)  +  cos  <}>x  (ex  +  ez  ny) | 


=  -90° 


in  equation 


(4.57) 


With  the  equations  for  a  gyroscope  and  the  equations  developed  in 
Chapter  II  for  a  rigid  body,  it  is  now  possible  to  obtain  the  equations  of  motion  for 
a  vehicle  containing  any  configuration  of  gyroscopes.  In  order  to  intelligently  choose 
a  configuration,  it  is  necessary  to  first  consider  the  nature  of  the  steady  state 
This  is  the  subject  matter  of  Chapter  V  . 


response. 
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chapter  v  -  steady  state  requirements  for  an  orbiting  vehicle 

Consider  a  vehicle  in  orbit  of  the  earth  that  contains  internally  mounted 
gyroscopes  for  damping  purposes.  If  no  disturbance  torques  are  applied,  steady 
state  will  eventually  be  attained.  If  the  vehicle's  orbit  is  circular,  and  the  gyroscope 
configuration  is  chosen  such  that  the  resultant  torque  exerted  by  the  gyroscopes  on 
the  body  at  steady  state  is  zero,  the  body  coordinate  system  will  become  aligned  with 
the  reference  coordinate  system;  the  body  will  then  have  an  angular  velocity  of  fl  j  . 
If  the  orbit  is  elliptical,  and  the  gyroscope  configuration  is  chosen  such  that  no  resul¬ 
tant  torque  is  exerted  on  the  body  at  steady  state,  0,0  and  their  derivatives  will 

7  x  z 

vanish  but  0^  and  its  derivatives  will  not  vanish  .  This  is  because  an  elliptical  orbit 
requires  the  vehicle  to  have  a  variable  angular  acceleration.  Since  the  restoring 
torques  are  produced  through  the  mechanism  of  differential  gravity  and  gyro^c  opic 
cross -coupling  of  the  vehicle' s  orbital  rate,  it  is  necessary  that  an  angular  deviation 
appear  about  the  axis.  It  should  be  noted  that  a  vehicle  in  an  elliptical  orbit  can 

never  attain  complete  alignment  with  the  reference  coordinate  system.  If  a  less 
restrictive  gyroscope  configuration  is  employed,  the  gyroscopes  will,  at  steady  state, 
exert  a  torque  on  the  vehicle.  In  order  to  counter  this  torque,  steady  state  will  be 
reached  with  the  body  axes  offset  from  the  reference  coordinate  system,  so  that 
counterbalancing  differential  gravity  and  gyroscopic  cross  coupling  torques  are 
produced.  Therefore,  in  steady  state,  the  best  that  can  be  expected  is  that 


ex*  ex  =  ez  =  ez  =  0  (5-D 

ex*  0  \  f  0  \  *  0  K* 0  K* 0  <5-2> 


Even  though,  for  an  arbitrary  gyroscope  configuration,  an  offset  will  occur  from  the 
reference  coordinate  system,  the  offset  can  be  determined  in  advance,  permitting 
the  orientation  of  data  gathering  equipment  to  be  adjusted  to  compensate  for  the 
angular  error.  It  is  therefore  anticipated  that  permitting  the  vehicle  to  have  a  small 
offset  from  the  reference  coordinate  system  will  not  produce  any  serious  problems. 
At  steady  state,  it  is  also  necessary  for  all  of  the  gyroscopes  in  the  configuration  to 
be  at  rest  with  respect  to  the  body;  that  is,  the  angular  velocity  and  acceleration  of 
each  gyroscope  gimbal  with  respect  to  the  body  be  zero.  Therefore, 


(5-3) 


for  each  gyroscope  in  the  configuration  at  steady  state.  The  desired  orientation  of 
each  gyroscope  at  steady  state  will,  of  course,  depend  on  the  vehicle  design.  The 
question  immediately  arises  as  to  the  nature  of  the  torque  that  must  be  supplied  to 
each  gyroscope  to  maintain  it  at  rest  in  its  desired  orientation  at  steady  state,  and 
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how  this  torque  will  be  supplied. 

A)  The  General  Steady  State  Torque  Equations 

When  a  gyroscope  is  mounted  within  an  orbiting  vehicle,  it  is  subjected  to 
an  angular  velocity  and  acceleration  due  to  the  vehicle1  s  orbital  rate-  If  the  gyro¬ 
scope  is  at  rest  with  respect  to  the  body,  the  angular  velocity  and  acceleration  of  the 
gyroscope  gimbal  with  respect  to  inertial  space,  is  the  same  as  the  angular  velocity 
and  acceleration  of  the  body  with  respect  to  inertial  space.  For  a  vehicle  with  an 
arbitrary  gyroscope  configuration  at  steady  state  in  an  elliptical  orbit,  the  angular 
velocity  and  acceleration  of  the  body  with  respect  to  inertial  space  can  be  obtained 
from  equations  (2.14)  through  (2.19)  by  setting 


e 

x 


0  = 
x 


e 

z 


»  n  =  n 


X  z 


=  0 


(5.4) 


t*)  i  /  =  CO  / 

b/l-S.  g/l.s. 
hJi-s.  g/i.s. 


9,  Vb^'VV  jb-  Qy  kb  <5'5> 

...  ,  1 

G  n  L  +  (0  )  j,  -  e  Q  k,  (5.6) 

zyo  y  y  Jb  x  y  b  1  ' 


The  components  of  the  torque  on  the  gyroscope,  necessary  to  produce  this  angular 
velocity  and  acceleration  at  steady  state,  can  be  obtained  along  the  gyroscope  axes 
from  equations  (4.24),  (4.25)  (4.26)  and  along  the  body  axes  from  equations  (4.29), 
(4.  30),  (4.  31)  by  setting 


G  =  0  -r  e  = 

X  X  Z 


(5.7) 


The  required  components  oftcrque  on  the  gyroscope  along  the  gyroscope  axes  are 

T  =  I  fsin  <J>  (0  +n  )  +  sin  A  cos  A  0  fi  (5.8) 

g  g  l  y  y'  yy  yz  x  y  9 

g  r 

■I*  (cos  <t>  cos  A  8  ft  1  +  H  [sin  A  cos  A  (0  +  ft  ), 

y  z  z  yj  r  {  x  rz  y  y 

+  (cos  A  cos  *  -  sin  A  sin  A  sin  d>  )  0  Q  -  (cos  A  sin  A 

x  ^y  x  _  Y  z  x  y  '  ^x  ^y 

+  sin  A  cos  A  sin  A  )  0  ft  1 
Tx  Yy  Yz'  z  yj 


T 


I  [cos  A  cos  A  (*0  +  ft  )  -  (cos  A  sin  A  sin  A 

g  t  x  z  y  y  x  y 

+  sin  A  cos  A  )  0  0  -  (cos  A  cos  A  sin  A 

Yx  Yy'  x  y  '  Yx  y  Yz 


sin  4>  sin  A  )  6 


I  f  -  sin  A  cos  A  (0  +  0  )  -  (cos  A  cos  A 

g  l  Tx  Yz  '  y  y  x  Y 

•  sin  A  sin  A  sin  A  )  ©  0  +  (cos  A  sin  A 

Yx  Yy  Yz'  x  y  '  Yx  Yy 

+  sin  <t>x  cos  <(>y  sin  <|>z)  6z  flyj  +Hr  [sin  <l>z  <0y  +  Oy) 

+  sin  4»y  cos  <t>z  ex  ny  ♦  cos  <{>y  cos  <)>z  ez  (yj 


(5-9) 


(5.10) 
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and  along  the  body  axes  are 

T  =  I  6  n  +  H  [(sin  6  cos  <J>  +  cos  6  sin  6  sin  d>  ) 

g:  g  z  Y  r  L'  Yx  Ty  Yx  Ty  Yz' 

(0  +  ft  )  +  cos  <J>  cos  (f>  0  fill 

y  y  Yx  Yz  x  yj 


(5.11) 


r  =1  (0  +  ft  )  +  H  [  (cos  4>  cos  <(>  sin  <j> 

g  g  y  y'  r  Yx  Ty  Tz 

Ju  *  • 

-  sin  4>  sin  <j>  )  0  ft  -  (sin  <)>  cos  <f>  +  cos  <)>  sin  <J> 

x  y  x  y  x  y  x  y 


(5.12) 


sin  <b  )  0 
z  z 


By] 


=  -  I  0  ft  +  H  fcos  d>  cos  <(»  G  f2  +  (cos  6  cos  d>  sin  d)  (5.13) 

Of  T  V  T»  7  V  '  ^  '  • 


g  x  y 


z  z  y 


-  sin  <b  sin  d>  )( 0  +0)1 

x  YyM  y  y'  J 

B)  Specialization  for  several  restrictive  cases. 

It  is  advantageous  to  specialize  these  equations  for  somewhat  more  restric¬ 
tive  conditions. 


1)  Consider  a  satellite  in  an  elliptical  orbit  with  a  gyroscope  configuration 
chosen  such  that  the  components  of  the  resultant  torque  exerted  by  the  gyroscopes  on 
the  body  along  the  and  axes  are  zero.  In  steady  state,  the  components  of 

torque  on  each  gyroscope  in  the  configuration  can  be  obtained  along  the  gyroscope 
axes  from  equations  (5.8)*  (5.9)*  (5.10)  and  along  the  body  axes  from  equations 
(5.11),  (5.12),  (5.13)  by  setting 

0x  =  0z  =  0  (5.14) 


The  required  components  of  torque  on  the  gyroscope  along  the  gyroscope  axes  are 

T  =  I  sin  <b  (0  +  ft  )  +  H  sin  6  cos  4>  (0  +  ft)  (5.15) 
g.  g  ^z  '  y  y'  r  Yx  Tz  y  y' 


T  -I  cos  (b  cos  d>  (0  +  0  )  (5.16) 

g;  g  *  z  '  y  y' 

Jg 

T  =  -  I  sin  <f>  cos  d)  (0  +  ft  )  +  H  sin  4>  (6  +  ft  )  (5. 17) 

gk  g  Yz  y  y'  r  rz  y  y' 

g 


and  along  the  body  axes  are 


H  (sin  A  cos  <b  +  cos  d> 
r  '  ^x  y  Tx 


•  (5.18) 

sin  <b  sin  <b  )  (0  +  ft  ) 

^y  Yz  y  y ' 


jb 


i  (e  +  )  (5.19) 

8  y  y 

Hr  ^COS  tx  COS  8in  ~  8in+x  8in<t,yX^r+ny)  <5‘2°) 
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2)  Consider  a  satellite  in  a  circular  orbit  with  a  completely  arbitrary  gyroscope 
configuration  employed.  In  steady  state  a  resultant  torque  will,  in  general,  be 
exerted  on  the  body  by  the  gyroscopes.  Therefore,  0^,  0^  and  0^  will,  in  general 
not  be  zero.  In  steady  state,  the  components  of  torque  on  each  gyroscope  in  the 
configuration  can  be  obtained  along  the  gyroscope  axes  from  equations  (5.8),  (5.9) 
(5.10)  and  along  the  body  axes  from  equations  (5.  11),  (5.  12),  (5  - 13)  by  setting 

0  =  0  =  O  =  0  (5.21) 

y  y  y 

The  components  of  the  torque  on  the  gyroscope  required  along  the  gyroscope  axes  are 

T  =  H  |  sin  A  cos  A  O  +  (cos  d>  cos  A  -  sin  A  sin  A  sin  A  )  0  ft 
g.  r  l  ~X  y  '  ~y  jC  ~y  2.  x  y 

lg  -I 

-  (cos  A  sin  A  +min  6  cos  A  sin  A  )  0  ft  T  (5. 22) 

^x  Yy  Tx  Yy  Yz'  z  yj  ' 


(5.23) 


T  =  H  f sin  A  cos  (|>en+  cos  A  cos  6  0  0+  sin  <b  ft  1  (5.  24) 

gR  r  [  y  x  y  Ty  Yz  z  y  yj 


and  along  the  body  axes  are 


r  =  H  r  (sin  A  cos  A  +  cos  A  sin  A  sin  A  )  ft 
g  r  l  Tx  >  Ty  Yz'  y 

*i_ 


+  cos  4>  cos  4>  0 


z  X 


nv] 


(5.25) 


(cos  A  cos  A  sin  A  -  sin  <J>  sin  A  )  0  ft 
l '  x  y  z  x  ^y  x  y 

(sin  A  co 8  A  +  cos  A  sin  A  sin  A  )  0  ft  V 
Yx  Ty  Yx  y  Yz'  z  y  j 


(5.26) 


T  =  H  f  cos  A  cos  A  0  ft  +  (cos  A  cos  A  sin  A 
gk  r  l  x  z  z  y  x  Yy  Tz 

b 

sin  A  sin  A  )  ft  1 
x  y  yj 


(5.27) 


3)  Consider  a  satellite  in  a  circular  orbit  with  a  gyroscope  configuration  chosen 
such  that  no  resultant  torque  is  exerted  by  the  gyroscopes  on  the  body  at  steady  state* 
The  equations  for  this  case  can  be  obtained  from  the  general  equations  (5.8)  through 
(5-  13)  by  setting 

0=0=050=0=0=0  (5.28) 

x  y  z  y  y  y 
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The  required  components  of  the  torque  on  the  gyroscope  along  the  gyroscope  axes  are 

T  =  H  sin  $  cos  <f>  O  (S  -  29 ) 

g.  r  Yx  Yz  y  '  ' 

g 

T  =0  (5.30) 

®  i 


T  =  H  sin 

e 

and  along  the  body  axes  are 

T  =  H  (sin 

\  r 


T  =0 

gi 

Jb 


<l>  n 

*  y 


<b  cos  <j>  +  cos  A  sin  A  sin  A  ) 

x  Yy  x  Ty  z  y 


(5.  31) 


(5.  32) 

(5.33) 


T  =  H  (cos  A  cos  A  sin  A 
g,  r  yx  Yy 

Kb 


sin  A  sin  A  )  U 
x  y  y 


(5.34) 


For  a  circular  orbit,  and  steady  state,  an  equivalent  form  for  the  torque 
on  each  gyroscope  in  a  configuration,  is,  of  course,  given  by 

Tg  =  «y  *  Hr  (5.  35) 

The  resultant  torque  on  all  the  gyroscopes  for  a  circular  orbit  at  steady  state,  is 
therefore  given  by 

E  (ny*Hr)  =  nyx  £ht  (5.36) 

The  resultant  torque  on  the  body  is  the  negative  of  the  resultant  torque  on  the  gyro  - 
scopes.  Therefore,  at  steady  state,  the  resultant  torque  on  the  body  will  be  zero, 
when  either  Y  =  0  (the  vectors  form  a  closed  polygon)  or  the  ^  H  *s  along 

n 

y 

At  steady  state,  it  is  necessary  for  the  required  to  que,  predicted  by  the 
equations,  to  be  present  on  each  of  the  gyroscopes  of  the  configuration.  Since  none 
of  the  gyroscopes  are  constrained  about  their  precession  axes,  it  will  be  necessary 
to  utilize  a  gyroscope  torquer  to  furnish  the  i^  component  of  the  torque  for  each 
gyroscope  for  which  this  component  is  not  zero.  For  a  circular  orbit,  the  required 
torque  will  be  a  constant.  For  an  elliptical  jrbit,  it  will  be  necessary  to  supply  a  time 
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varying  torque  which  will  undoubtedly  require  the  use  of  programmed  digital  equip¬ 
ment-  It  would  be  interesting  to  consider  what  would  happen  if  gyroscope  torquers- 
are  not  utilized.  For  a  circular  orbit,  the  precession  angle  of  each  gyroscope  in 
the  configuration,  and  the  body  angular  deviations,  will  adjust  themselves  so  that  the 
resultant  torque  on  the  body  will  be  zero.  If  desired,  the  values  of  the  precession 
angles  and  body  deviation  angles  can  be  obtained  by  first  setting  equation  (5.22)  equal 
to  zero  to  obtain  a  relation  between  the  precession  angle  and  the  body  deviation  angles, 
and  then  equating  the  negative  of  the  resultant  torque  applied  to  the  gyroscopes  (the 
resultant  torque  is  obtained  by  summing  equations  (5.25),  (5.26),  (5.27)  for  each 
gyroscope  in  the  configuration  )  to  the  torque  on  the  body,  equations  (2.  52),  (2.  54), 
(2.56)  after  setting  all  the  derivatives  equal  to  zero. 


-St 

8H 

-  Yj  T 

gj 

Jt 

• 2  V 


-  i  )  4  e  n  c 

yb  x  y 


i  )  e  u 

zb  y  y 


-i  )  n  n 

xb  z  y 


(5*37) 

(5-  38) 

(S-  39) 


If  the  orbit  is  elliptical,  the  gyroscope  will  never  reach  a  state  of  rest  with  respect 
to  the  body.  Therefore,  energy  will  be  constantly  dissipated  due  to  gyroscope 
damping.  The  only  sources  of  this  energy  are  the  kinetic  and  potential  energy  of  the 
satellite.  Invariably,  as  time  progresses,  the  orbital  parameters  will  therefore  be 
alte  red 

C)  Choice  of  a  Gyroscope  Configuration 

Consider  briefly  the  problem  of  choosing  a  gyroscope  configural ion  for  a 
vehicle  in  a  circular  orbit-  The  configuration  must,  of  course,  be  chosen  to  provide 
damping  about  all  vehicle  axes.  (This  aspect  of  the  problem  will  be  considered  in 
detail  in  Chapter  VI.  )  If  possible,  the  gyroscope  configuration  should  be  chosen  so 
that  no  resultant  torque  is  exerted  on  the  body  due  to  the  gyroscopes  at  steady  state- 
If  this  cannot  be  achieved,  it  will  be  necessary  to  determine  the  steady  state  offset 
angles  from  equations  (5.  37),  (5.  38),  and  (5-  39)-  (The  precession  angle  of  each 
gyroscope  is  specified  by  the  configuration.  )  The  torque  which  the  torque r  must 
supply  to  each  gyroscope,  is  given  by  the  i^  component  of  the  torque  on  the 
gyroscope. 

T  =  T 

torque r  g. 


(5.40) 
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It  is  also  necessary  for  each  gyroscope  in  a  configuration  to  be  in  a  state  of  stable 
equilibrium.  The  nature  of  the  equilibrium  can  be  easily  established  for  each 
gyrost  ope  by  substituting  the  value  of  4>x  specified  by  the  design  into  the  derivative 
of  equation  (5-  22)  with  respect  to  4>x  . 


(5  41) 


u  r 

Tg.  =  Hr  [cos  ’t’x  co8  ^  ny  -  (sin<f>xC08  4>y +cos  <bx  sin<by  sin<|>z)0xny 
‘g 

-  (cos  4>  cos  <(>  sin<J>  -  sin  4>  sin<f>  )  0  fi  1 
x  y  z  x  y '  z  yj 


If  equation  (5  41)  is  positive,  a  state  of  stable  equilibrium  exists.  If  equation  (5.41) 
is  negative,  a  state  of  unstable  equilibrium  exists*  The  explanation  for  this  test  is 
as  follows:  Equation  (5.41)  is  the  derivative  of  the  torque  required  to  maintain  the 
gyroscope  at  any  specified  position  at  steady  state  If  the  derivative  is  positive,  at 
the  specified  position,  a  positive  angular  displacement  of  the  gyroscope  will  require 
a  larger  torque  than  is  being  supplied  to  maintain  the  new  position.  A  negative 
angular  displacement  will  require  a  smaller  torque  than  is  being  supplied  to  maintain 
the  displacement.  Hence,  in  either  case,  a  resultant  torque  will  be  established  such 
that  the  gyroscope  will  return  to  its  original  position-  If  the  derivative  .s  negative^, 
a  positive  angular  displa<  nment  will  require  a  smaller  torque  than  is  being  supplied, 
and  a  negative  angular  displacement  will  require  a  larger  torque  than  is  being 
supplied.  In  either  case,  the  resultant  torque  will  be  established  such  that  the 
gyroscope  will  continue  to  move  in  the  direction  displaced. 


In  order  to  obtain  a  better  understanding  of  the  principles  involved,  an 
example  will  be  considered.  Consider  the  configuration  shown  in  Fig  4-*  where  a 
single  gyroscope  is  mounted  with  its  precession  axis  along  the  axis  of  the  vehicle. 

The  satellite  design  requires  the  spin  axis  of  the  gyroscope  to  be  at  an  angle  <j>x  with 
the  Y,  axis.  It  will  be  assumed  for  convenience  that  the  satellite  is  in  a  circular 

D 

orbit.  It  is  immediately  apparent  that  if  $  is  not  zero,  steady  state  will  be  estab¬ 
lished  with  the  body  offset  from  the  reference  coordinate  system.  If  the  vehicle  was 
aligned  with  the  reference  coordinate  system  at  steady  state,  the  angular  momentum 
vector  ,  will  no*  be  along  0  The  angular  offset,  at  steady  state,  can  be 

obtained  by  setting  4>  =  4>  =  0  in  equations  (5.25),  (5. 26)  and  (5  27)  and  substituting 

y  z 

in  equations  (5.  37),  (5-  38)  and  (5.  39) 


-  H  (sin  *  n  +  cob  <b  0  0  )  =  -  (I  -  I  )40x 


(5.42) 


r’  Tx  y  rx  x  r  '  zb  yb 


x  y 


Hr  9in  +x  e„  °y  =  3(Ixb  *  I*b)  6y  °y 


2 


(5.43) 
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•  H  cos  <J> 
r  x 


0  ft 

z  y 


=  (I  - 1  ) 

yb  xb 


e  n 

z  y 


(S.44) 


Solution  of  these  equations, 


0  = 
x 


yields  the  offset  angles  of  the  vehicle  at  steady  state. 
H  sin  <i> 


r  x 


-I  )  4fi  -  cos  d>  H 
xb'  y  Yx  r 


(5.45) 


If 

H  cos  6  -  (I  -  I  )  *  0  (5.46) 

r  x  xb  yb  y 

then 


0=0=0 

y  * 


(5  47) 


If  the  left  side  of  equation  (5.46)  is  equal  to  zero,  0^  and  0^  are  indeterminate- 
The  torque  which  must  be  supplied  by  the  gyroscope  torquer  is  given  by 


T  =  H  n  (cos  ()>  0  +  sin  *  )  (5.48) 

torquer  r  y  x  x 

If  a  gyroscope  torquer  is  not  utilized,  steady  state  will  be  established  with  4>x  =  0 
and  the  vehicle  aligned  with  the  reference  coordinate  axes. 

If  another  gyroscope  is  added  to  the  configuration  of  Fig.  4-5  so  that  the 
configuration  of  Fig  6~2,  with  4>z  =  0  results,  steady  state  will  be  established  with 
the  vehicle  axes  aligned  with  the  reference  coordinate  axes,  since  for  this 
configuration,  the  resultant  angular  momentum  vector,  at  steady  state,  is  along 
fty  .  This  configuration  is  known  as  a  "  VM  configuration  and  will  be  discussed 
further  in  Chapter  VI. 


In  Chapter  VI,  the  results  of  Chapters  II,  IV  and  V  will  be  combined  to 
obtain  the  equations  describing  the  dynamics  of  an  orbiting  vehicle  with  internally 
mounted  gyroscopes.  Several  configurations  will  be  considered  that  provide 
damping  about  all  of  the  vehicle  axes.  For  simplicity,  the  treatment  will  be 
restricted  to  vehicles  in  a  circular  orbit. 
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CHAPTER  VI  -  THE  DYNAMICS  OF  AN  ORBITING  VEHICLE  WITH  INTERNALLY 


MOUNTED  GYROSCOPES  FOR  A  CIRCULAR  ORBIT 

The  dynamics  of  an  orbiting  vehicle  in  a  circular  orbit  with  internally 
mounted  gyroscopes  will  be  developed  for  several  gyroscope  configurations.  There 
is  no  contention  that  the  configurations  chosen  for  presentation  are,  in  any  sense, 
optimal;  they  are  presented  to  illustrate  the  principles  involved.  The  gyroscope 
configurations  chosen  for  analysis,  will  provide  indirect  damping  about  all  three 
of  the  vehicle  axes.  It  was  noted  in  Chapter  III  that  a  smaller  settling  time  could  be 
obtained  through  the  utilization  of  indirect  damping  than  with  direct  damping.  It  is 
for  this  reason  that  all  the  configurations  presented  are  chosen  to  provide  indirect 
damping.  In  fact,  if  direct  damping  was  desireable,  it  would  be  merely  necessary 
to  construct  the  vehicle  in  two  parts,  coupled  viscously,  since  the  mechanism  of 
gyroscope  precession  is  not  involved.  It  has  been  shown  in  the  literature,  and  should 
be  expected  from  the  results  of  Chapter  III  of  this  report,  that  direct  vehicle  damping 
will  result  in  excessive  settling  time. 

The  dynamics  of  an  orbiting  vehicle  with  internally  mounted  gyroscopes 
can  be  developed  by  combining  the  results  obtained  in  Chapter  II  for  an  orbiting  rigid 
body  and  the  results  obtained  in  Chapter  IV  for  a  gyroscope  under  various  constraints. 
The  torque  exerted  on  the  body  is  equal  to  the  vector  sum  of  the  external  disturbance 
torques  (exclusive  of  the  differential  gravity  torque)  exerted  on  the  body  and  the 
torques  exerted  by  the  gyroscopes  on  the  body.  For  a  symmetrically  designed  gyro¬ 
scope,  there  is  no  differential  gravity  torque  on  the  gyroscope,  hence  the  torque 
exerted  by  the  gyroscope  on  the  body  is  the  negative  of  the  torque  exerted  on  the 
gyroscope 


where  Tft  is  the  external  torque  applied  to  the  body  less  the  differential 

gravity  torque. 

T^*  is  the  total  torque  on  the  body  less  the  differential  gravity  torque. 

T  is  the  torque  on  the  n*1  gyroscope  of  the  configuration. 

®n 

A)  Analysis  of  a  Configuration  Employing  Gyroscopes  with  Precession  axes 

along  and  . 

Fig.  6^1  shows  a  configuration  employing  two  gyroscopes  which  will 
provide  damping  about  all  vehicle  axes.  One  of  the  gyroscopes  is  oriented  with 
its  precession  axis  along  the  Z,  axis  (precession  angle  denoted  by  <|>  )  and  the  other 

D  X, 
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Fig.  6-1  Configuration  Employing  Gyroscopes  with  Precession 
Axes  Along  and 
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with  its  precession  axis  along  the  X,  axis,  (precession  angle  denoted  by  4>  ).  If 

x2 

gyroscope  torquers  are  not  utilized  for  both  gyroscopes,  at  steady  state,  the  body 
will  be  aligned  with  the  reference  coordinate  axes  and  the  steady  state  values  of 

4>  and  A  ,  A  and  4>  ,  will  be  zero.  For  such  a  condition,  indirect 

1  X2  1  x2 
o  o 

damping  will  not  be  attainable  about  the  axis.  To  obtain  indirect  damping  about 
the  Y^  axis,  it  is  necessary  that  at  least  one  of  the  steady  state  gyroscope  precession 

angles,  <{>  or  4>  ,  be  non-zero.  In  general,  for  arbitrary  values  of  4>  and 

1  X2  X1 
°  °  o 

4>x  indirect  damping  about  the  axis  is  provided  by  the  Z^  oriented  gyroscope, 

o 

indirect  damping  about  the  Yfe  axis  is  provided  by  both  the  Xfe  and  oriented 

gyroscopes,  and  indirect  damping  is  provided  about  the  Zfc  axis  by  the  Xfe  oriented 

gyroscope.  For  arbitrary  values  of  <|>  and  <)>  ,  a  gyroscope  torque r  will  be 

X1  x2 
o  o 

required  for  each  gyroscope. 

The  angular  deviations  of  the  vehicle  as  a  function  of  <J>  and  <J>  at 

X1  x2 
o  o 

steady  state  can  be  calculated  from  equations  (5.25),  (5.  26),  (5.  27),  (5.37),  (5.  38) 
and  (5.  39).  Solution  of  these  equations  yields: 


H  sin  A 
r 


o  4(1  -  I  )  n  -  H  (cos  A  4  cos  A  ) 

zb  yb  y  r  X1  x2 


(6.2) 


-H  sin  A 
r  x 
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0  = 
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o  (I  -  I  )  ft  +  H  (cos  <|>  +  cos  <}>  ) 

yb  ^  y  r  X1  x2 
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(6.3) 


0  = 


H  (sin  $  0  sin  A  0  ) 

r  z  x.  x  ' 

2  o  1  ° 

_ o _ o 

3(1  -1)0 

xb  V  y 


(6.4) 


The  torque  which  each  gyroscope  torquer  must  supply  can  be  determined  from 
equation  (5.40).  For  the  X^  oriented  gyroscope,  the  torque  required  is  obtained  from 
equation  (5.22)  with  $  =  ^  *  0  . 

T 

torquer 


=  H  fi  (Bin  4>  +  co.  *  6 i  )  (6.5) 

r  y  x2  *2  *, 
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For  the  oriented  gyroscope*  the  torque  required  is  obtained  from  equation  (5.22) 

with  4>  =-90°  and  A  =  0* 

'  y  Yz 


torque r 


=  H  ft  (sin  <J>  +  cos  $  0  ) 

if  ’  TV  Tv  T 


r  y 


x,  z 
1  o 
o 


(6.6) 


whe  re 


0^  and  0z  are  obtained  from  equations  (6.2)  and  (6.3) 


o  o 

The  values  of  4>  and  <J>  must,  of  course,  be  chosen  so  that  the 

X1  x2 

o  o 

gyroscopes,  at  steady  state,  are  in  a  state  of  stable  equilibrium.  To  check  the  nature 
of  the  equilibrium  it  is  merely  necessary  to  substitute  the  steady  state  deviation 
angles  (equations  (6.2),  (6.  3),  (6.4))  into  equation  (5.41)  for  each  gyroscope  and 
examine  the  sign. 


For  stability 
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The  dynamic  performance  of  the  vehicle  with  internally  mounted  gyroscopes 
can  be  obtained  from  the  body  equations  for  a  circular  orbit  (2.52),  (2.54),  (2- 56), 
from  the  equations  of  an  oriented  gyroscope  (4.46),  (4.47),  (4.48),  from  the 
equations  of  a  Z^  oriented  gyroscope  (4.54),  (4.55),  (4.56)  and  from  equation  (6.1). 
The  resulting  equations  are: 
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=  (i  +2i)e  +3(1  -i  )  e  n  + 1 
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g  y 


y  y 
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♦x2  (ez  •  exV  (MO) 


-  4>  (0  +  0  fi  )  I  -  H  sin  <|>  (0  -  0  n  )  +  sin  4>  (0  +  <>  n  ) 

x^x  z  y  i  r  XjZ  x  y  x^  x  z  y 


sin  <(>  +  <j>  sin 


X1  x2 


r  =  (i  +21  )(e  -  e  n  )  +  (i  -i  )(en  +00) 

z,_  g'  z  X  y'  yh  x^'  x  y  z  y  ' 


z  y 


+  I 


4>  -  <L  (®v  +ft„)  +H  A  cos  A 

i  2  y  y  J  [x2 


+  8 in  <(>  (8  +  )  +  (cos  A  +  cos  <)>  )  ( 0  +00) 


X1  y 


X 2'  'X  Z  y' 


(611) 


The  equation  of  constraint  for  the  X,  oriented  gyroscope*  relating  <(>  to  the  body 

D  x 2 

deviation  angles,  can  be  obtained  from  equation  (4.49),  where  Tc  is  equal  to 
T 

torquer^  given  by  equation  (6.5) 


torquer 


=  D(*>  +  I  (0  +  0  0  +  A  )  +  H  lsin  4>  (6+0) 

x2  g  x  z  y  x2  r|_  x2  y  y 


-  cos  A  (0-00) 
tx2  z  x  y 


(6.12) 


The  equation  of  constraint  for  the  oriented  gyroscope  can  be  obtained  from 

equation  (4.  57)  where  T  is  equal  to  T  given  by  equation  (6  6) 

rq  r^ 


torquer 


-  D<i>  +  I  (b  -  e  n  +  %  )  +  H  [sin  <b  (0  +  U  ) 

7.  xi  *  z  x  y  xi  rL  xi  y  y 

V] 


+  cos  A  (9+0 
x^  x  z 


(613) 


The  equations  developed  in  this  section  characterize  the  steady  state  and 
dynamic  performance  of  an  orbiting  vehicle  for  the  configuration  of  Fig.  6-1  for 

arbitrary  4>  and  <(>  .  Since  the  equations  are  nonlinear*  the  dynamic  performance 

X1  x2 
o  o 

of  a  vehicle  containing  this  configuration  can  best  be  analyzed  using  an  analog  or 

digital  simulation.  It  will  be  necessary  to  choose  values  of  4>  ,  4>  and  find  a 

X1  x2 
o  o 

suitable  vehicle  geometry  so  that  both  adequate  settling  time  is  attained  and  the 
vehicle  angular  deviations  at  steady  state  are  small.  Indeed,  for  large  angular 
deviations,  the  approximations  made  in  Chapter  11  will  probably  invalidate  the 
re  stilts. 
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Fig.  6-2  "V"  Configuration  in  the  Plane 
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B)  Analysis  of  a  "V"  configuration  in  the  Y^  plane.  . 

Fig.  6^2  shows  two  gyroscopes  with  their  precession  axes  in  the 
plane  forming  a  "V"  configuration.  This  configuration  will  provide  indirect  damping 
about  all  three  of  the  vehicle  axes.  At  steady  state,  one  of  the  gyroscopes  will  have 
its  spin  axis  above  the  Y^  plane  by  an  angle  .  This  gyroscope  will,  for 

convenience,  be  referred  to  as  "gyro  I"  .  The  othei^gy  rose  ope  will,  at  steady  state, 
have  its  spin  axis  below  the  Y^  plane  by  an  angle  <t>^  .  This  gyroscope  will 

*o 

for  convenience,  be  referred  to  as  "gyro  II"  .  The  precession  angle  of  gyro  I  will 

be  denoted  by  4^  and  the  precession  angle  of  gyro  II  will  be  denoted  by  4>  •  At 

1  x, 

steady  state 


<t>  =  4>  (6.14) 

xi  xi 


4>  =  360°-*  (6.15) 

x2  X1 

o 

The  angular  deviations  of  the  vehicle  as  a  function  of  4>  at  steady  state 

xi 

can  be  calculated  from  equations  (5. 25),  (5.26),  (5.27),  (5.37),  (5 .  38)  and  (5 .  39). 

H  cos  4>  cos  4>  6  *  2(1  -I  )  6  n  (6.16) 

r  x,  z  x  z,  y,  '  x  y 

1  o  b  yb  o  7 

o 


2  H  cos  <i>  sin  4>  ® 

r  Tx.  x 

Xo  ° 


-  I 

xb 


n 

y 


(6.17) 


2  H  (cos  4>  cos  4>  6  +  cos  4>  sin  4>  ) 

r  x.  z  z  x«  z 

1  o  1 

o  o 

=  (I  -1)0  n 
*b  vb  *0  y 


(6.18) 


Solution  of  these  equations  yields  the  offset  angles  of  the  vehicle  at  steady  state. 


2  H  cos  4>  sin  4> 
r 

e  = - 2 - 

zo  (I  -  1  )  n  •  +  2H  cos  4>  cos 


(6.19) 


yb  ^  y 


‘i 


if 


n  cos  *  cos  <)>  -  2(1  -  I  )  n  £  0 

r  *1  2  2b  V  y 

O 


(6.20) 


then 


(6.21) 
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If  the  left  side  of  equation  (6.20)  is  equal  to  zero,  0  and  0  are  indeterminate. 

xo  yo 

Note,  that  if  4>  =  0,  0  will  also  be  zero  and  hence,  at  steady  state,  the  vehicle 

z  z 

o 

will  be  perfectly  aligned  with  the  reference  coordinate  axes. 

The  torque  which  each  gyroscope  torquer  must  supply  can  be  determined 
from  equation  (5.40).  For  gyro  I,  the  torque  required  is  obtained  from  equation 


(5 . 22)  with  <jy  =  0 


(6.22) 


torquer^  r  y 
r 


=  H  O  (sin  A  cos  <j>  +  cos  <b  0  -  sin  A  sin  A  0  ) 

1*  tY  if  tY  _  7,  7 


X,  X 

1  o 

o 


For  gyro  II,  the  torque  required  is  obtained  from  equation  (5.22)  with  Ay  =  0  and 

4>  =  360°  -  <b 


1 


=  H  fi  (cos  A  0  -sin  A  cos  A  +  sin  A  sin  A  0  ) 


(6.23) 


torquer^  r  yx  Tx,  x_  '  Tx,  Tz  Tx,  Tz  z_ 


1  o 
o 


'1 


Since  0X  =  0 
o 


torquer^ 


=  H  Q  sin  A  (cos  A  -  sin  4>  0  ) 

r  T  r  y  x«  z  z  z 


(6.24) 


H  fi  sin  A  (sin  A  0  -  cos  A  ) 


torquejjj  r  y  ^ 


(6.25) 


where  0z  is  obtained  from  equation  (6.19). 
o 

The  value  of  A  must,  of  course,  be  chosen  so  that  both  gyroscopes, 

xi 

o 

at  steady  state,  are  in  a  state  of  stable  equilibrium.  To  check  the  nature  of  the 
equilibrium^  is  merely  necessary  to  substitute  the  steady  state  deviation  angles 
(6. 19)  and  (6.21)  into  equation  (5.41)  for  each  gyroscope  and  examine  the  sign.  For 
stability  of  both  gyroscopes 


cos  A  cos  A  + 

X1 


2H  cos^<b  sin^A 
r  TXj  z 

o 


(I  -  I  )n  +  2H  cos  A  cos  A 
Yb  xb  *  r  X1 


>0  (6.26) 


If  the  precession  axes  of  the  gyroscopes  are  along  X^(AZ  =  0)>  the  condition  for 
stability  reduces  to 


cos  A  >  0 

Tx 

X1 

o 


(6. 27) 
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Clearly, the  gyroscopes  will  be  in  a  state  of  stable  equilibrium  when 

-  90° <<j>  <  +  90°  (6.28) 

X1 

o 

The  dynamic  performance  of  the  vehicle  with  internally  mounted  gyroscopes 
can  be  obtained  from  the  body  equations  for  a  circular  orbit  (2. 52),  (2.  54),  (2. 56), 
from  the  equations  for  a  gyroscope  whose  precession  axis  is  in  the  plane 

(4.  34),  (4.  35),  (4.  36)  and  from  equation  (6.1).  The  resulting  equations  are: 


=  (I  +  21  )  <e  +  e  oj  +  (i  -  i  )  <e  n  -  40  n  ) 


z  y 


z  y 


x  y 


(6.29) 


+  I 
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+  H 
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z  X1  x2 


z  z  x  y 


(0  +  O  )  (sin  A  +  sin  Aj  -cos  A  (0  -  0  O  )  (cos  A  +  eo«  6  ) 

y  y  z  «  x  y  xj  x2 


+  sin  A,  (A„  sin  +  A„  8in  A„  ) 

z  Xj  Aj  x  ^  x  2 


T  =  (I  +  21  )  0  +  3(1  -1)00+1 

y,  g  y  xk  *K  y  y  8 

Jb  b 


+  A  )  (6-  30) 

x2 


+  COS 

+  sin 


[8in  ^Z  \ 

A,  (©,  -  0  O  )  (A  +  A,  )  -  H  cos  A.  (A,  sin  A,  +  A_  sinA_  ) 

2  2  x  y  Xj  x2  J  r,L  X1  X1  x2  X2 
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%  zb  *  z  x  y  Vb  xb  X  y  z  y 


(6.  31) 


+  I 
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+  H 


sin  A.  (©,  +  ©,  n„)  (A,  +  A  )  -  COS  A  (©  +  O  )(A  +  A  ) 

1  x2  2  y  y  X1  x2  J 


1 


2  '  X  z  y  ■x1  'x2'  rz’  y  y”TXj  x, 
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The  equation  of  constraint  for  gyro  I,  relating  its  precession  angle  to  the  body 
deviation  angles,  can  be  obtained  from  equation  (4.  37),  where  Tc  is  equal  to 

Ttorquer  given  in  «9uation  (6.  24) 

T 

torque r 


a  +  I  U 

ri  xi  g  L 

+  H  [sin  4> 
r  1  xi 


0  sin  4>  +  cos  4>  (0  +  0  fi)+4> 

y  ~z  z  x  z  y  Tx 
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(6.32) 


cos  4>  (0  +fl  )  -  cos  4>  (0-00) 
Yz  '  y  y'  Yx1  z  x  y' 
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~ Z  X  Z 
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The  equation  of  constraint  for  gyro  II  is  identical  with  that  for  gyro  I  (equation  (6.  32)) 

with  <j>  .  replaced  by  6  and  Tx  replaced  by  T 

Xj  r  7  torquer^  r  7  torquer^. 

Although  the  development'  for  the  "V"  configuration  in  this  section  was 
for  equal  and  opposite  displacement  angles  from  the  plane,  the  development 

could  have  been  easily  extended  to  the  case  where  the  displacement  angles  are 
arbitrary.  In  fact,  an  extension  to  a  "V"  configuration,  where  ‘he  precession  axes 
lie  on  any  line  relative  to  the  body  coordinate  system,  can  be  easily  accomplished 
from  the  equations  developed  in  Chapters  II,  IV  and  V  in  a  manner  similar  to  that 
used  in  this  section. 

C)  A  Configuration  Employing  a  Single  Gyroscope  to  Provide  Three  Axis 

Vehicle  Damping  * 

Fig-  6-3  shows  a  configuration  employing  a  single,  one  degree  of  freedom, 
gyroscope  whose  precession  axis  is  located  in  the  plane.  This  configuration 

will  provide  damping  about  all  three  of  the  vehicle  axes.  Since  only  one  gyroscope 
is  used,  no  subscripting  is  necessary  therefore,  the  gyroscope  precession  angle  will 
be  denoted  merely  by  <t>^,  and  the  value  of  the  gyroscope  precession  angle  at  steady 
state  will  be  denoted  by  <)> 

o 

The  angular  deviations  of  the  vehicle  as  a  function  of  <j>x  at  steady  state 
can  be  calculated  from  equations  (5.25),  (5.26)  and  (5.27)  with  =  0  and  from 
equations  (5. 37),  (5.  38)  and  (5.  39).  Solution  of  these  equations  yields: 

H  sin  <b  cos  <J> 
r  x  y 

e  a - 2 -  (6.  33) 

xo  4(1  -  I  )f2  -  H  cos  4> 

*b  yb  y  r  xo 


H  sin  <b  sin  6 
r  x  v 

_ o _ 1 

(I  -  I  )Q  +  H  cos  <j> 

yb  V  y  r  5 


H  sin  6  (sin  d>  6  +  cos 

r  x  y  x 

n  '  n 


3(i  -  i  )n 

xb  zb  y 


*>  e  ) 

y  z 
'  o 


(6.34) 


(6.35) 
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Fig*  6-3  Configuration  Employing  a  Simple  Gyroscope 
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The  torque  which  the  gyroscope  torque r  must  supply  can  be  obtained  from 

equation  (5.22)  with  <)>  =  0  and  from  equation  (5.40) 

2  (6.36) 

T,  =  H  n  (sin  A  +  cos  A  cos  A  0  -  cos  A  sin  A  0  ) 

torque  r  r  y  Yx  Yx  y  x  x  Yy  z 

n  7  o  o7o  o7o 

where  G  and  0  are  obtained  from  equations  (6.  33)  and  (6.  34) 
x  z 

o  o 


The  value  of  4>  must,  of  course,  be  chosen  so  that  the  gyroscope,  at 
o 

steady  state,  is  in  a  state  of  stable  equilibrium.  To  check  the  nature  of  the  equilibrium 
it  is  merely  necessary  to  substitute  the  steady  state  deviation  angles  (equations 
(6-  33),  (6.34),  (6.35)  into  equation  (5.41)  and  examine  the  sign.  For  stability 


cos  A  cos  A  -  sin  d>  cos  A  0  +  sin  A  sin  A  0  >  0  (6.  37) 

Yx  Yz  Yx  Yy  x 1  Yx  Yy  z  ' 

o  o  7  o  o  7  o 

The  dynamic  performance  of  the  vehicle  can  be  obtained  from  the  body 
equations  for  a  circular  orbit  (2.52),  (2.54),  (2.56),  from  the  equations  for  a  gyroscope 
constrained  to  rotate  in  the  plane  (4.  38),  (4.  39),  (4.40)  and  equation  (6.1)  . 

The  resulting  equations  are: 
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The  equation  of  contraint  relating  the  gyroscope  precession  angle  to  body 
deviation  angles  can  be  obtained  from  equation  (4.41)  where  Tc  is  equal  to 

Ttorquer  given  by  e9uation  36) 


T.  =  D<fr  +  I 

torque r  g 


cos  6  (6  + 

Yy  x 


0  ft  )  -  sin  <)>  (0  - 

z  y  y  z 


e  n  )  + 

x  y' 


+  H 


sin  <|>  (0  +  )  -  cos  6  cos  d>  (0  -  0  fl  ) 

X  V  V  v  T»r  7  v  w' 


-  cos  6  sin  (6  +  0 

Tx  y  x  z 


v] 


x  y' 


(6.41) 


The  procedure  involved  for  developing  the  equations  describing  the  statics 
and  dynamics  of  an  orbiting  vehicle  for  an  arbitrary  configuration  should  be  clear 
from  the  examples  that  have  been  presented.  The  equations  developed  in  Chapters 
II,  IV  and  V  can  be  used  to  analyze  any  configuration  in  a  circular  orbit.  Further¬ 
more,  many  of  the  equations  can  be  applied  to  a  vehicle  in  an  elliptical  orbit 
Chapter  VII  will  summarize  the  results  thus  far  obtained  and  will  indicate  the  areas 
in  which  further  study  is  warranted. 
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Chapter  VII  -  CONCLUSION 

The  analyses  given  in  this  report  are  valid  if  the  following  assumptions 
are  satisfied: 

1.  Attitude  deviations  are  small  so  that  the  equations  of  motion  for  the  body  can 
be  linearized. 

2.  With  the  exception  of  the  internally  mounted  gyroscopes,  the  satellite  comprises 
a  rigid  body  with  no  moving  internal  parts. 

3.  Control  is  desired  about  the  principal  axes  of  the  vehicle. 

4.  The  earth  is  assumed  to  be  a  homogeneous  sphere,  and  therefore  the  gravitational 
field  is  inverse  square. 

5.  All  gyroscopes  are  assumed  symmetric. 

It  is  further  assumed  in  the  latter  part  of  Chapter  V  and  Chapter  VI  that  the  orbit  is 
circular. 

With  the  above  assumptions  the  following  has  been  accomplished: 

1.  The  equations  of  motion  of  a  rigid  body  in  an  elliptical  orbit  were  obtained  with 
the  inclusion  of  the  gravity  gradient  torque. 

2.  An  analysis  of  a  simple  configuration  employing  gyroscope  damping  was 
considered  to  furnish  an  insight  into  the  mechanism  of  gyroscope  damping  and  its 
prope  rties. 

3.  The  equations  of  a  symmetrical  one  degree  of  freedom  gyroscope  mounted  in  an 
arbitrary  position  within  the  vehicle  in  an  elliptical  orbit  were  obtained  and  several 
specialized  cases  were  cited. 

4.  The  steady  state  requirements  for  an  orbiting  vehicle  containing  interna ilv 
mounted  gyroscopes  were  considered-  In  the  latter  part  of  this  analysis  the  vel.ir  le' <■ 
orbit  was  assumed  circular. 

5.  The  equations  of  motion  for  a  vehicle  containing  internally  mounted  gyroscopes 
were  obtained  for  several  gyroscope  configurations  for  a  circular  orbit. 

There  are  several  areas  of  investigation  that  have  been  neglected  in  this 
work  and  should  be  the  topic  of  future  study  in  this  field. 

1.  The  effects  of  large  angular  deviations  should  be  considered  since  at  the  time 
of  insertion  of  the  vehicle  into  its  orbit  large  angular  deviations  will  invariably  occur. 
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2.  The  effects  of  an  non-spherical,  non-homogeneoufl  earth  should  be  considered. 

3.  The  effects  of  moving  internal  components  should  be  investigated.  The 
acceleration  and  deceleration  of  internally  moving  parts  (such  as  tape  recorders) 
will  produce  reaction  torques  on  the  vehicle.  Motion  of  components  may  also  cause 
the  moments  of  inertia  and  location  of  the  vehiclers  principal  axes  to  become  time 
varying. 

4.  The  analysis  should  be  extended  to  include  non- symmetrically  designed 
gyroscopes. 

5.  A  computer  study  of  several  gyroscope  configurations  is  essential  to  obtain 
both  an  optimum  vehicle  geometry  and  optimum  gyroscope  configuration. 
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APPENDIX  A 

Derivation  of  Euler1  a  Equations 


For  any  system  of  particles,  the  angular  momentum  of  that  system  about 
any  point  in  space  (see  Figure  A-"l) 


Fig.  A-l  System  of  Particles 


is  given  by 


where 


n 

h  =  X  r.  x  m.  v. 
i=i  1  11 

th 

m.  is  the  mass  of  the  i  particle 


(A.  1) 


r.  is  the  radius  vector  from  the  point  about  which  the  angular  momentum 

*  ^ 

is  to  be  obtained  to  the  l  particle. 

v.  is  the  velocity  of  the  i**1  particle  relative  to  the  point  about  which 
the  angular  momentum  of  the  system  is  to  be  obtained. 


This  definition  is  meaningful  regardless  of  whether  the  reference  point  be  fixed  or 
translating  with  respect  to  inertial  space* 


The  case  where  the  system  of  particles  comprises  a  rigid  body,  and  the  point 
about  which  the  angular  momentum  is  to  be  obtained  is  the  mass  center  of  the  body,  is 
of  particular  interest.  For  this  case,  the  velocity  of  an  arbitrary  (i**1)  particle  is 
given  by 


v.  “  w  x  r. 
i  i 


(A.  2) 
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where  <o  is  the  angular  velocity  of  the  rigid  body  with  respect  to  inertial  space 
is  the  radius  vector  from  the  mass  center  to  the  i^  particle. 

Substituting  into  the  general  equation  for  the  angular  momentum,  (A.l)  one  obtains  for 
the  angular  momentum  about  the  mass  center 

n 

hc.B.  =  ^  r  x  m.  («  x  r  )  (A.  3) 

®  i=l 

Using  the  vector  identitiy  for  a  triple  cross  product 

A  x  (B  x  C)  =  B(A.  C)  -  C(A.  B) 


one  obtains 


r.  x  m.  (to  x  r.)  =*  m. 
11  l  l 


wr.  -  r.  (u  •  r.) 


Resolving  r.  and  w  into  a  system  of  rectangular  coordinates 

r.  =  xi  +  yj  +  zk 

w  =  w.i+w.i+w.  k 
i  J  J  k 

and  performing  the  indicated  operations,  one  obtains 

h  =h.  i  +  h.j+h1  k 

c.g.  1  J  k 


where 


h,  =  A  to.  -  H  co  .  -  G  co. 
i  i  j  k 

h.=  -  H  co .  +  B  w.  •  F  w, 

J  i  J  k 


(A  - 4) 
(A. 5) 

(A. 6) 
(A.  7) 

(Art) 

(A. 9) 
(A.  10) 
(A.  11) 


h.  =  -  Gw.  -  Fw.+Cw. 

*  l  j  k 

Equations  (A.  9),  (A.  10)  and  (A.  11)  can  be  expressed  in  matrix  notation  by  representing 
the  angular  momentum  vector  and  the  angular  velocity  vector  as  column  matrices 

/  A  -  H 


h.  / 

1 

h. 

J 

= 

!  h  i 

\  w  \ 

-  H 


B 


-  G  -  F 


\ 

/.«.  \ 

1 

CO  . 

J 

W 

(A.  12) 


The  first  term  on  the  right  side  of  the  equation  is  referred  to  as  the  inertia  tensor. 
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Where 

A  = 

Z  m.  (y2  +  z2) 

(A.  13) 

B  = 

Z  mi  (*2  +  *2) 

(A. 14) 

* 

C  = 

Z  mj  <*2  +  y2) 

(A- IS) 

and 


are  the  polar  moments  of 


inertia  about  the 

i,  j  and  k  axes  respectively 

X 

u 

m.  xy 

(A.  16) 

O 

ii 

m.  xz 

i 

(A.  17) 

T  »  Z 

m.  yz 

(A. 18) 

the  i-j, 

.i-k,  and 

are  the  products  of  inertia  in 
j-k  planes  respectively. 

If  i,  j  and  k  are  body  fixed  vectors,  the  components  of  the  inertia  tensor  are  <.  onstanl*- 


It  is  shown  in  any  text  of  classical  mechanics  that  the  torque  exerted  on  a  body 
about  its  center  of  gravity  is  equal  to  the  rate  of  change  of  its  angular  moment  urn  about 
its  center  of  gravity  as  seen  from  inertial  space. 


where 


T  =  h„  „  (A.  19) 

eg.  c.g. 

6  h 

h  =  - +  wxh  (A  20) 

c*g*  6t  C#g# 


whe  re 


^/fit  denotes  a  partial  differentiation  in  which  the 
i-j-k  vectors  are  held  fixed, 


(*>  is  the  angular  velocity  of  the  i-j-k  coordinate 

system  with  respect  to  inertial  space  (whic  n  need  not 
not  be  that  of  the  body). 


It  can  also  be  shown  in  advanced  books  of  classical  mechanics  that  there  exists 
at  least  one  set  of  axes  for  which  the  product  of  inertia  terms  (H,  G,  F)  vanish  Such 
axes  are  known  as  principal  axes.  For  the  principal  axes,  the  inertia  tensor  is 
diagonal. 


II 


0 

0 


0 

B 

0 


(A.  21) 
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It  will  be  assumed  that  the  i-j-k  vectors  are  fixed  in  the  body  and  coincide  with  the 
principal  axes . 

Therefore 

h  =  Aco.i+Bco.j+Cco^k  (A- 22) 

eg.  i  J-  k. 


and 


•  #  •  ® 

n  =  A(o,  i  +  Bw.  j  +  Cu,  k 

Cg-  1  JJ  k 

+  cox  lAco^i  +  B(o^j+Cw^k| 


(A.  23) 


Expanding  and  equating  (A.  23)  component  wise  to  the  torque  yields  the  following 
equations : 


Ti  - 

A  co .  +  (C  -  B)  co  .  co, 
i  J  K 

(A.  24) 

Tj  - 

BcOj+(A-C)ioi<oc^ 

(A.  25) 

Tk  • 

$ 

C  <ov  +  (B  -  A)  co .  co  . 

K  1  J 

(A.  26) 

These  are  Euler*  s  equations. 
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APPENDIX  B 

Derivation  of  Gyroscope  Equations. 

A  gyroscope  consists  of  two  basic  components.  An  inertia  wheel  or  rotor 
which  rotates  at  a  high  speed  and  a  gimbal  to  which  the  rotor  is  mounted  and  which 
permits  angular  motion  of  the  rotor  with  respect  to  inertial  space.  (Refer  to  Fig*B-l). 

The  torque  on  the  rotor  and  gimbal  which  comprise  the  gyroscope,  about  the 
center  of  gravity  of  the  gyroscope  is  given  by  the  time  derivative  of  its  angular  momen¬ 
tum  about  its  center  of  gravity  with  respect  to  inertial  space- 


T 

g 


c* 


g* 


h 

gc.  g. 


(B.l) 


The  angular  momentum  of  the  gyroscope  can  be  considered  to  be  composed  oi  two 
parts.  One  due  to  the  angular  velocity  of  the  rotor  and  the  other  due  to  the  angular 
velocity  of  the  gimbal  .  The  component  of  angular  momentum  due  to  the  angular  velo¬ 
city  of  the  gimbal/  is  computed  by  considering  the  rotor  to  be  stationar/. 


*c.  g. 


Now 


5  A  r.  x  m.  v.  +  ).r. 

u  i  11.  u  i 

I  rotor 

),  r.  x  m.  v.  I 

u  i  i  i| 


x  m.  v. 
l  ii 


I  rotor 
H 


=  Hrjg 


J  W 
r  r 


(B  d) 

gimbal  phis 
station*  ry  roto r 


(B-  3) 


(B.4) 


where 


represents  the  polar  moment  of  inertia  of  the  rotor  about  its  axis 
of  rotation. 


u>  is  the  angular  velocity  of  the  rotor  about  its  axis  of  rotation- 


Er.  x  m.  v. 
l  l  i 


...  ,  =  A  u>  i+Bu.j+Cw, 

gimbal  plus  gig  g  S  g  k„ 

g  g  g 


(B5) 


|  stationary  rotor 
As  was  done  in  deriving  Euler’s  equations  in  Appendix  A 


*  hg 


eg. 


6h 

a  r£  +  w,  xh 


g  g 


(B.6) 


whe  re 


<*>=<*>.i+a>.  j+u.  k  (B-7) 

g  ig  8  Jg  %  \  « 

is  the  angular  velocity  of  the  gimbal  axes  with  respect 
to  inertial  space- 


Fig*  B-l  Single  Degree  of  Freedom  Gyroscope 
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Therefore 


d  ,  *  •  *  ♦ 

TT-h  s  A  u,  i  +  (H  +  B  w .  )  j  +  C  (*>. 

«  gc.g.  «  ’g  g  *  g  jg  '8  g  kgfg 


(B.8) 


♦  .  f 

w3  L 


(C  -  B 
8  g 


+  k 


)  co  *  co.  -  H  a).  +  j  (j 

i  k  r  k 

g  g  g  J  ^  L 

(B^  -  A  )  co.  co  .  +H  co .  u 
g  g  i  j  r  x  • 

L  *  g  Jg  gj 


(Ag  •  cg)  toi  wk  , 
8  8  g  g  J 


Grouping  terms  and  equating  component  wise  to  the  torque  one  obtains 

T  =*  A  co,  +  (C  -  B  )  co  .  co.  -  H  co. 
g.  gig  g  J„  k„  r  kc 


g  g 


g 


T  aH+Bu.  +(A-C)u.  w. 

%  r  8  Je  88  ’«  \ 

=C  co.  +(B  -  A  )  co .  co.  +H  <o, 

rr  n  Lr  '  n  rr  ’  «  <  «•  . 


g 


g 


g  k  '  '  g  g'  i  wj  ’  r  wi 
8  g  8  8  g  Jg  g 


(B.9) 

(B.10) 

(B.  11) 


Examination  of  these  equations  reveals  that  the  torque  on  the  gyroscope  is 
composed  of  two  parts.  The  first  part  is  given  by 


T  aA  i,  +(C  -B  )w,  w, 

%  *  ‘8  *  *  jg  kg 

T  *  B  u.  +  (A  -C)w.  u. 

%  8  J8  8  8  '8  kg 

T  =  C  <*>.  +(B  -A  )«.  co. 

\  *  kg  *  8  ■«  >t 


and  the  second  part  is  given  by 


T  =  -  H  «. 

g.  r  k 

n 


T  =  H 

gj  r 
Jg 


T  a  H  «. 

*k  r  1 
g 


(B.  12) 
(B.  13) 
(B. 14) 

(B.  15) 
(B.  16) 
(B.  17) 


g 


The  first  part  is  seen  to  be  a  restatement  of  Euler's  Equations,  (A,  24),  (A.  25),  (A.  26). 
It  describes  the  rigid  body  dynamics  that  do  not  depend  upon  the  angular  momentum 
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of  the  rotor.  The  second  part  of  the  equation,  which  is  known  as  the  gyroscopic 
torque  ,  is  provided  by  the  angular  momentum  of  the  rotor  .  If  the  angular  velocity  of 
the  rotor  is  assumed  to  be  constant  (Hr  »  0)  ,  the  equation  of  gyroscopic  torque  takes  a 
particularly  simple  form.  Namely 

T  =  wxH  (B.  18) 

g  r 

where  u>  is  the  angular  velocity  of  the  gimbal  with  respect  to  inertial 
space. 

Hy  is  the  angular  momentum  of  the  rotor. 

This  form  of  the  equation  is  extensively  used  in  the  literature. 


114 


APPENDIX  C 


Derivation  of  the  Gtavity  Gradient  Torque. 

Consider  a  rigid  system  of  particles  in  orbit  around  the  earth.  (See  Fig.  C*l) 


The  force  on  the  i**1  particle  is  given  by 


W.  =  -  m.  g.  — 
—  i  1  i  p. 


(C.  1) 


Where 


Now 


m.  is  the  mass  of  the  i**1  particle 


.th 


g.  is  the  acceleration  due  to  gravity  at  the  position  of  the  i 
particle . 

p,  is  the  vector  displacement  from  the  center  of  the  earth  to  the 

•th  *  *  | 
i  particle. 


Re 

1  ^ 


(C.2) 


where  g8  is  the  acceleration  due  to  gravity  at  the  surface  of  the 
earth. 

is  the  mean  radius  of  the  earth. 


Substituting  eq.  (C-2)intoeq.  (C.l)  yields: 


W.  =  -  m.  g  -  p. 

—  1  1  °8  3 


(C.3) 


The  torque  exerted  on  the  system  of  particles  about  its  center  of  gravity  due  to  the 
earth' 8  gravitational  field*  is  given  by  n 

T_dg=£iiX-*  <C-4> 


1 


where  is  the  gravity  gradient  torque  on  the  body  (differential  gravity 

—  torque) 

r.  is  the  vector  displacement  from  the  center  of  gravity  to  the 

.th  . 

i  particle* 


Substituting  eq.  (C.  3)  into  eq.  (C*4)  yields: 

l  y  m. 

^dg  =  -  8s  Re  -3 


(C5) 
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From  Fig.(Cwl)  it  is  seen  that 


p.  =  p  k  +  r. 

i-i  ro  —  —i 


(C- 6) 


where 


Now 


p  is  the  distance  from  the  center  of  the  earth  to  the  mass  center  of 
ro 

the  system  of  particles. 

k  is  a  unit  vector  outward  along  the  geocentric  vertical- 

3 


pi  3  =  (pi  •  £*)  "  7 


(C-7) 


Substituting  eq.  (C.6)  and  simplifying  yields: 


-3  -3 

pi  =  po 


2k.  r. 

r.  \2 

1  +  7-"1  + 1 

JL\ 

P 

ro 

» 

,  3 

7 


Since 


x-«i 

o 


Pi 


-3 


1  + 


2k-  r. 


3 

7 


Expanding  in  a  power  series  and  taking  only  the  first  terms  yields; 


-3  ^  -3 

Pi  r?  P0 


1  - 


3k  • 


(C.8) 


(C.9) 


(C-  10) 


(C.  11) 


Substituting  eq.  (C.6)  and  (C- 11)  into  eq.  (C.5)  and  expanding  the  vector  product  yields 

g  R  2  n  3g  R  2  n  (C- 12) 

Jig  =  ~TT~  -  X  i?j  mi-i  +  ~  r~  45j  mi<£ixi>Lri;y 


But 


T,  m.  r.  =  0 
i-I  »-* 


(C- 13) 


defines  the  center  of  mass  of  the  system.  Therefore 
T.  -38'R'2 


dg 3  -  s  ~  xi>  <ii  i> 


(C- 14) 
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For  a  circular  orbit 


P  n 
Ko  y 


g  R 
*8  6 


(C .  15) 


where 


Therefore 


gQ  is  the  acceleration  due  to  gravity  at  the  mass  center. 


g8  R  2 

— — =  n  c 


(C.  16) 


Substituting  eq.  (C.16)  into  equation  (C*14)  yields  for  a  circular  orbit. 


n 

T  .  =  3  n  2  Y.  m.  (r.  x  k)  (r.-k) 

_ dg  y  i=i  1  -1  -  -i - 


(C.  17) 


Let  X^»  Y^  and  be  the  coordinate  system  formed  by  the  body  principal  axes, 
with  corresponding  unit  vectors  i^,  and  .  Let  X,  Y  and  Z  be  the  reference 
coordinate  system  with  corresponding  unit  vectors  i,  j  and  k  .  (Refer  to  Chapter  II 
flections  A  and  B.  )  The  set  of  Eulerian  angles  relating  the  body  coordinate  system 
to  the  reference  coordinate  system  is  denoted  by  0^,  0^  and  0^  .  The  transformation 
matrix  relating  the  two  coordinate  systems  is  given  in  equation  (2.5)  .  (Refer  to 
Chapter  II,  section  C.) 


Expressing  the  vector  i\  in  terms  of  body  coordinate  8  yields : 

ii =  *b  H> +  yb  jb  +  zb  kb 

Using  the  transformation  matrix  (2.5)  yields: 

k  =  -  0  i,  +  0  j,  +  k 
—  y  b  x  Jb  b 

Substituting  equations  (C.19)  and  (C.18)  into  eq.  (C.14)  yields: 

2  „ 


T  = 
Adg 


3  g„  R 


l 

i=l 


m. 

l 


(C.18) 


(C.  19) 


(C. 20) 


S>  (Vb  '  ex  zb  >  +  ^  ('  *b  '  ey  zb) 

+  kb(exXb  +  6y  V]  [-Vb  +  exyb+,b] 
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i 


Expanding  eq.  (C.20)  and  noting  that  the  product  of  inertia  terms  vanish  (see 
eqe .  (A.  16),  (A.  17),  (A.  18)  yields : 


3*S  R 


n1  r 

-  & m*  I1 


H,  W-  V  ♦  Jb  6v '\2- %2| 


4  kb  9,  6y  <%*  - 


(C-  21) 


The  component  will  be  neglected  since  it  involves  the  product  of  two  small 
angles.  Substituting  eqs.  (A. 13),  (A.  14)  and  (A.  15)  into  eq.  (C.  21)  yields  the  final 
result  for  the  gravity  gradient  torque. 


R 

_ _ e 

3 


2 


)+jb 


e  (i 
y  a 


) 


(C-  22) 


Substituting  eq.  (C.16)  into  eq.  (C*  22),  one  obtains  the  corresponding  equation  for  a 
crircular  orbit. 


3  n 


Led 

D  X  Z 


•V4i» 


e  (i 


(C. 23) 


I 
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